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CHAPTER 1 
INTRODUCTION 


1 . 1 Background 

The ability to measure certain ocean wave characteristics 
over large areas and long periods of time could conceivably have 
significant Impact upon open ocean and coastal activities. Ship 
routing, search and rescue operations, meteorological research and 
recreational activities are just a few of the areas where quick and 
reliiible sea state information is desired. 

With the advent of satellite altimetry, it is now possible 
to estimate many ocean wave characteristics with a degree of accu^ 
racy which equals or exceeds previous techniques. Specifically, the 
altimeter on board the Geodynamics Experimental Ocean Satellite 
(GEOS-3) sampled the radar return waveform from which the ocean sig- 
nificant waveheight (SWH) could be inferred. The ocean significant 
waveheight. Which is sometimes referred to as H-1/3, is defined as 
the average of the one-third highest waves in a long sequence of 
waves observed at a point (Neuman and Pierson, 1966). 

Satellite significant waveheight measurements taken in the 
more remote regions of the earth's oceans are of particular interest 
since local estimates of significant waveheight are made almost 
entirely from ships at infrequent intervals. The small number of 
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ships passing through such remote regions, combined with the 
inherent inaccuracies of human ”eyeball” estimates, make global sig- 
nificant waveheight measurements from satellites highly desirable. 

The GEOS-3 estimates of significant waveheight were made 
on a near-global basis and in near-real-time. Although preliminary 
estimates of the measuremer'; accuracy have been made (McMillan and 
Roy, 1977f and Parsons, 1977 and July, 1979), a formal determination 
of th&' accuracy of these estimates has yet to be presented. 

1 . 2 General Description of GEOS -3 Spacecraft 

The GEOS-3 spacecraft (see Figure 1.1) was launched from 
the Air Force Western Test Range on April 9, 1975, as part of the 
National Geodetic Satellite Program, with the specific objectives of 
improving man's knowledge of the earth's gravitational field, the 
size and shape of the terrestrial geoid, the deep ocean tides, sea 
state, current structure, crustal structure, solid earth dynamics 
and remote-sensing technology (Stanley, 1979). The spacecraft 
orbited the earth with a period of 101.8 minutes in a near-circular 
orbit with an inclination to the equator of 115**. Data were col- 
lected and archived from the satellite until November 30, 1978, when 
active operations were terminated, 

GEOS-3 telemetry data were acquired from four sources: 
the NASA STDN (Space Tracking and Data Network) VHF telemetry sites. 
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FIGURE 1,1, QEOS-3 SPACECRAFT CUTAWAY 
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FIGUKE 




TOTAL 

NUMBER 

OPERATIONAL DATES 

STATION 

OF GOOD 
PASSES 

START 

STOP 

Merritt lilandi, Florida (MIL) 

list 

April 10, 1975 

Dec. 1, 1978 

Rotman, North Carolina (ROS) 

641 

April 10, 1975 

Dec. 1, 1978 

Winkfiald, United Kingdom (WNK) 

1221 

April 10, 1975 

Dec, 1, 1978 

Bermude Itlandt (BOA) 

398 

April 10, 1975 

Dec, 1, 1678 

Madrid, Spain (MAO) 

122 

April 10, 1975 

Dec, 1, 1978 

Aiceniion (ACN) 

457 

April 10, 1975 

Dec, 1, 1978 

Johanneiburg (BUR) 

118 

April 10, 1975 

Oct . 31, 1975 

Guam (GWM) 

410 

April 10, 1975 

Dec. 1, 1978 

Orroral, Australia (ORR) 

482 

April 10, 1976 

Doc, 1, 1978 

Hawaii (HAW) 

706 

April 10, 1975 

Doc, 1, 1978 

Fairbanks, Alaska (ULA) 

1051 

April 10, 1975 

Dec, 1, 1978 

Gti^^ltone, California (GOS) 

486 

April 10, 1975 

Doc. 1, 1978 

Quito, Ecuador (QUI) 

383 

April 10, 1975 

Dec, 1, 1978 

Santiago, Chile (AGO) 

395 

April 10, 1975 

Dec, 1, 1978 

Tananarive, Madagascar (TAN) 

12 

April 10, 1975 

July 11, 1975 

Mahe, Seychelles (MAH) 

79 

March 1, 1976 

May 9, 1976 

Rosman (via ATS-6 94° West) 

78 

April 10, 1975 

June 12, 1975 

Rosman (via ATS-6 140° West) 

347 

Sept. 3„ 1976 

Dec. 1, 1978 

Madrid, Spain (via ATS-6 34° East) 

34 

May 26, 1975 

Oct, 22, 1976 


TABLE 1,1, NASA TELEMETRY STATIONS 
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OHiUli 1 . . .. 

OF POOR QUALITY 


TOTAL 

NUMSEn OPERATIONAL DATES 

OF GOOD 


STATION 

PAT,SES 

START 

STOP 

Htmdon, Virginia (HER| 


April 20, 1676 

July 1, 1976 

Parth, Adftralia iAUSi 

82 

July 14, 1076 

Aug, 17, 1976 

Tafuna, Samoa (TAF| 

120 

Aug, 12, 1976 

Nov, 30, 1976 

Shamva, Ala^a |SHM) 

234 

Aug. 18, 1676 

Oac. 6, 1676 

Napiar, Ntw Zaaland (NEZ| 

8B 

Aug, 26, 1676 

Oct, 24, 1676 

Eaatar liland |EAS| 

222 

Nov. 8, 1676 

Fab, 1, 1676 

Ba)a, Maaico |BAJ) 

137 

Oac. 12, 1676 

Fab, 1, 1976 

Kargualan laland IKEG| 

137 

March 2, 1676 

May 16, 1676 

Cmoi lilanda ICOCI 

113 

March 13, 1676 

May 28, 1676 

Falkland Itiarrdi <FLK) 

186 

April 21, 1676 

June 20, 1676 

Canary Itlandt |CY|| 

128 

July 4, 1676 

Sept, 7, 1676 

Natal, Brazil <NAT| 

B3 

July 14, 1676 

Nov, 16, 1676 

Triitan Da Cunha (TOK) 

32 

Aug. 3, 1976 

Oct. 14, 1976 

Caribou, Malna (CAR I 

640 

Sapt. 6, 1676 

Dac, 11, 11)76 

Kwaialein (KWA) 

37 

Jan. 16, 1677 

April 17, 1677 

Saattia, Vi/aihington (SEAl 

166 

Jan, 18, 1677 

April 3, 1677 

Rangiroa, Tahiti (TAH> 

73 

Jan, 22, 1677 

Fab. 23, 1977 

Towntviilt, Auatralia (TOW| 

72 

March 16, 1677 

July 20, 1677 

Cyprui (CYP) 

161 

April 18, 1677 

July 22, 1677 

Maha, Saychallat (MAH) 

101 

Aug, 10, 1877 

Nov, 13, 1977 

Salalah, Oman (SAL| 

233 

Aug, 16, 1677 

Nov, 26, 1677 

Pratoria, South Africa (PRE| 

206 

Aug. 27, 1977 

Nov, 30, 1677 

Adak Island, AlaAa (ADAI 

B9 

Jan. 20, 1978 

March 26, 1678 

Kourou, French Guiana (KOU) 

62 

May 5, 1678 

Juna 20, 1678 

Papeeta, Tahiti (TAHI 

B 

Dac, 16, 1677 

Jan. 20, 1678 

Pinang, Malaysia (MALI 

32 

Mays, 1678 

July 2, 1678 

Okinawa (OKI) 

48 

May 2. 1977 

July 24, 1977 


TABLE 1.2. DEFENSE MAPPING AGENCY TELEMETRY STATIONS 



the portable DOD (Depar tment of Dofenac) VHF racilltles, the NASA 
STDN S-band telemetry altea and the NASA ATS-6 aatellite telemetry 
data relay link (sea Tables 1.1 and 1.2 and Figure 1,2, reproduced 
from Stanley and Dwyer, 19.80). 

The primary instrument on board the GEOS-3 spacecraft was 
a radar altimeter, developed for NASA by the General Electric Cor- 
poration. This altimeter operated at a single frequency of 13.9 GHz 
and transmitted 100 radar pulses per second. The pulses transmitted 
by the altimeter were reflected from the earth's surface and 
received by the spacecraft. The time interval required for a pulse 
to make the round-trip could be used to determine the altitude of 
the spacecraft above the mean earth surface. In addition, the slope 
of the return pulse received at the spacecraft could be used to 
determine the characteristics of the surface. 

The GEOS-3 altimeter was Instrumented with 16 sample and 
hold gates (see Section 1.3), which provided information about the 
shape and amplitude of the return waveform. This information could 
be used to determine a number of interesting and useful parameters, 
including the spacecraft attitude, water/land and water/ice boun- 
daries, surface wind speed and significant waveheight. Significant 
waveheight was determined through analysis of the return waveform at 
the NASA Wallops Flight Center (WFC), Wallops Island, Virginia, dur- 
ing the entire active mission of the spacecraft. 
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The GEOS-3 spacecraft structure, which was patterned after 
the GEOS-2 structure, was an octahedron topped by a truncated 
pyramid. Extending from the pyramid was a gravity-gradient boom 
with an end mass, which was extended after insertion into orbit. 
The spacecraft was oriented in a stable gravity-gradient attitude, 
where the direction was defined by the gravitational force acting on 
the spacecraft (approximately radial). 

The radar altimeter was mounted on the octrahedron oppo- 
site from the pyramid and the boom (facing radially inward toward 
the earth^s surface). The altimeter antenna was required to be 
aligned within 1.2 degrees of the local vertical at all times. This 
pointing precision was maintained by the gravity-gradient boom and 
end mass configuration and by a constant speed, angular momentum 
wheel. Estimates of the pointing angle error have been made using 
the radar return waveform and were found to comply with the 1.2 
degree pointing requirement (McMillan, November I960). 

Also included in the Instrument package were coherent and 
non-coherent C-Band transponders, laser retroflectors, doppler 
transmitters and S-Band instrumentation for earth tracking and 
satell ite-to-satell ite exper iments . 

1 . 3 Lineage of the GEOS -3 Altimeter 

The Skylab S-193 altimeter was the first in the series of 
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satellite altimeters that were planiKHi to pr oyr essl voly achieve the 
goal of 10 oin revSolutiorr iir the iiotellite altitude above the ocean 
surface. That experimental aUi.ieter was designed primarily to 
obtain the radar measur ements trecessary for designing improved 
altimeters. The UK0S"“d altimeter^ seeorul in the series of satellite 
altimeters, was the first altimeter system applied to global opera- 
tion. The Advanced Applications Flight Experiments (AAFE) Altlrrre- 
ter, an aircraft system which first collected data in October 1975, 
Was a developmental effort directed at bridging the technology gap 
between the capabilities of the (iEOS-3 altimeter and the rather 
stringent requirements imposed on the Seasat altimeter, as well as 
providing surface truth in support of the Seasat altimeter calibra- 
tion activities. The Seasat spanonraffc, launched June 26, 1978, 
carried the third in the series of satellite altimeters and 
represented the first attempt to achieve 10 cm altimetric precision 
from orbit. It was conceptually identical bo the AAFE altimeter. 
The Seasat spacecraft also carried other sensors which were dedi- 
cated to ocean applications (IEEE Journal of Ocean ic Engineering, 
April I960). 

P esc ription of the flEOS-.S Altimotiy;;. 

'Die GEOS-3 altimeter operated in two distinct data gather- 
ing modes; namely, the intensive or short-pulse mode and the global 
or long-pulse mode. The global mode was the original mode designed 





for GEOS. However, the Natiotiai Oeoanlc and AtmoaphorAo Administra- 
tion (NOAA' and the Department of Defense (DOD) also decided to sup- 
port the intensive mode, which used pulse compression (at that time 
unproven over oceans). The intensive mode proved to be so accurate 
and reliable that it was used for the vast majority of the GEOS-3 
data segments. 

The global mode transmitter consisted of a magnetron tube 
with a 200 nanosecond pulse width capable of measuring height to a 
precision of \ meter. It also provided a measurement of the back- 
scatter signal intensity but did not use the 16 sample and hold 
gates that the intensive mode employed. Tl^orefore, it provided 
minimal information aboyt the return pulse shape. 

The Intensive mode transmitter consisted of a traveling 
wave tube ('IVT) with u pulse eompresslou to 12 nanoseconds and pro- 
vided improv/ed accuracy derived from a shorter pulse width, a more 
stable transmitter and waveform samples. Tlie receiver was coherent 
with both transmitters. 

The intensive mode gate timing and positioning is illus- 
trated in Figure 1.3 (reproduced from Mofmeistor and Keeney, 1977). 
As can be seen from that figure* four tracking gates wore positioned 
at various points along the waveform in order to define the waveform 
shape. These tracking gates were: 
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FIGURE 1-3- tMTEMtlVE llOOe GATE TIIIWIG 
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The Noise Gate, This gate had a width of 200 ns with a leading 
edge MOO na in front of the leading edge of the return pulse 
and provided a reference for measuring the amplitude of the 
return waveform, 

2, The Ramp Gate. Tiiis gate had a width of 12,5 ns and was 

located on the leading edge of the return pulse. It was used 
as a time reference for the other gates, 

3, The Plateau Gate, This gate had a width of 12,5 ns and was 

located 62.5 ns after the leading edge of the return pulse. It 
was used to determine the magnitude of AGO attenuation neces- 
sary to normal i*e the waveform. 

M. The Attitude/Specular Gate, iTiis gate had a width of 200 ns 

with a leading edge 700 ns after the leading edge of the 

Waveform. It was used to provide information about antenna 
bandwidth effects, wind speed and pointing angle. 

In addition, a waveform sampling system of sample and hold 
gates was included in the intensive mode to provide detailed infor- 
mation concerning the sha^«/ of the altimeter return waveform. Six- 
teen sample and hold gates were provided and positioned fixed in 
time with respect to the tracking gates (see Figure 1.3). The width 
of the gates was designed such that the leading edge of the average 
impulse response of the ocean surface could be obtained for 



waveheights vip to 10 meters to witnin * 20 %. As will be demon- 
strated in this investigation, this accuracy re<4uirement was met and 
exceeded. More specific information conerniag the function and 
positioning of the gates can be obtained from Hofmeister, et i«l 
(1976). 


The mean return waveforms were normalized using an 

automatic gain control (AGO system. While tracking, the AGC 
attenuation was adjusted to hold the average plateau gate output to 
a constant value. The other tracking gatea and the sample and hold 
gates were adjusted by the same amount. The return waveform could 
then be analyzed to determine some of the characteristics of the 
reflecting surface (the sea surface). 

Telemetry data werci transmitted in three formats for each 
of the two altimeter modes. The telemetry (TM) formats are summar- 
ized below: 

TH #1 - 2.0 seconds per data record with altitude data fre- 

quency of 10 per second 

TM #2 - 3.2 seconds per data record with altitude data fre- 

quency of 10 per second 

TM #3 - 3.2 seconds per data record with altitude data fre- 

quency of 100 per second 



,lt should be noted that, due to the power constraints of 
the spaoecraft, the altimeter was not operational at all times. 
Altimeter data segments were generally limited to tracks over water 
and were scheduled several weeks in advance of the actual data 
acquisition. This complexity of scheduling was not merely an incon- 
venience, since data from local weather disturbances, such as hurri- 
canes, could not be acquired unless the altimeter had been scheduled 
weeks earlier to take data at that time and place. 

For more information concerning the GEOS-3 spacecraft 
hardware subsystems see GEOS -C Hisaion Plan , NASA, 197*1. 

I. 5 GEOS -3 Significant Waveheight Estimation 

The initial GEOS-3 significant waveheight estimation algo- 
rithm was developed by G. S. Hayne (Hayne, 1977) and programmed by 

J. D. McMillan (McMillan, 1975) at WFC for use as a quality control 
check on the GEOS-3 altimeter preprocessing software. Various 
improvements and refinements of the algorithm by McMillan and Roy 
(1977) led to the achievement of a high degree of agreement between 
the GEOS-3 estimates and ship-based measurements of SWH. 

Eventually, the SWH estimate was distributed to GEOS-3 
principal investigators as an integral part of the altimeter data 
set. Various experiments performed at WFC by McMillan and Roy 
(1977) and by Fedor, et al (1979) indicated that the estimate of 



significant wavehelght produced by the GKuS-3 altimeter preprocesa-- 
ing software compared favorably with other estimates and measure- 
ments of sea state. 

The first comprehensive comparison of GEOS-3 significant 
waveheight estimates with independently derived sea state Informa- 
tion was presented by McMillan and Roy (1977) at the GEOS-3 Final 
Investigators' Meeting in New Orleans, Louisiana, In November, 1977. 
That study presented several variations of the original significant 
waveheight algorithm and included results obtained using different 
convergence criteria, different risetime coefficients and different 
ARS timing and amplitude biases (see Chapter 2). In that investiga- 
tion, the accuracy of the GEOS-3 SWH estimate was detei mined to be 
55 cm. 


Tlie significant waveheight estimate proved to be so useful 
that, in 1978, NASA established the GEOS-3 Near-Real-Time Data Sys- 
tem (McMillan, 1978) for disseminating significant waveheight and 
wind speed estimates. This system employed the significant 
Waveheight algorithm which was developed by McMillan and Roy. 

In the near-real-time system, the GEOS-3 altimeter data 
were acquired through the NASA Space Tracking and Data Network loca- 
tions (see Table 1.1) or the ATS-6 satellite and transmitted to the 
NASA, Goddard Space Flight Center in Greenbelt, Maryland, where they 
were buffered to magnetic tape in real-time. The buffered data were 



then transmitted to the Computer Sciences Corporation INFONET center 
in Beltsville, Maryland, where significant waveheight was computed 
and made available to user-supplied terminals on a call-up basis, 

During February, 1976, near-real-time GEOS-3 significant 
waveheight estimates In the North Atlantic Ocean were closely moni- 
tored and compared to other significant waveheight measurements pro- 
vided by the Spaceflight Meteorology Group of the National Weather 
Service and the Navy Fleet Numerical Weather Central. In Parsons 
( 1976 ), comparisons were made between the data sets, and It v<as 
found that the inherent consistency of the GEOS-3 data makes the satellite 
product the best representation of the true sea state. 

Even though the significant waveheight estimation algo- 
rithm as developed by McMillan and Roy was employed on thousands of 
GEOS-3 passes, no definitive study has been undertaken to establish 
the accuracy of the estimates. This Investigation will establish 
the accuracy of the significant waveheight estimates contained on 
the GEOS-3 data tapes (which are archived at the National Oceano- 
graphic Data Center in Suitland, Maryland) and will present a global 
atlas of all of the significant v;aveiieight data processed during the 
GEOS-3 mission. 

1.6 Applications of Satellite Significant Waveheight Estimates 

The global coverage and established accuracy of the altlm- 
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eter signiricant waveheight measurementa ^nake the data set prefer'- 
able to other sources Of information, which rely almost entirely 
upon infrequent and sometimes unreliable ’’eyeball” measurements that 
are typically in error by two to three meters. Consequently, the 
formulation of a global atlas of SWH has potential use by government 
and private industry to augment the information currently being 
implemented in climatic models. These models typically describe the 
oceans and the atmosphere in three dimensions but require Informa- 
tion concerning the air/sea interface. While the near-real-time 
data discussed previously would be extremely useful in short-term 
forecasts, the historical data presented in a global atlas could 
help define more accurately the relationship between the ocean and 
atmospheric models for which data have already been taken. 

For example, the U. S. Navy and the Coast Guard have 
expressed interest in the significant waveheight data to assist in 
the planning and execution of search and rescue operations. Again, 
the near-real-time data have been requested for real-time search 
operations and will be available from future spacecraft. The his- 
torical contour information, hewever, is also highly desirable for 
the study of the statistical probability of locating an object at 
sea in a given area at a given time of year. 

The Bureau of Land Management, as well as the petroleum 
industry, has expressed an interest in the data for use in the 
design of offshore facilities. Of particular interest is the 
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relationship between significant waveheight and the structural 
stress endured by offshore drilling equipment. An historical atlas 
of this type could help determine the statistical probability of 
high stress situations and, therefore, the design of equipment for 
particular areas. For example, drilling equipment in the North Sea 
has completely different structural requirements than does drilling 
equipment in the Baltimore Canyon (see, for example, McMillan, 
December, 1980). 

According to NOAA, significant waveheight data in some 
parts of the world are virtually nonexistent. This problem is par- 
ticularly acute in the southern hemisphere, where the lack of data 
has prevented the establishment of well-defined shipping lanes. 
Tlirough the use of a significant waveheight atlas, the accuracy of 
waveheight measurements in the southern hemisphere would be greatly 
enhanced. The result would be a much more accurate data base from 
which to determine shipping lanes on a seasonal basis. 

The U. S. Coast Guard is charged with the responsibility 
of controlling oil spills in the ocean areas near the United States 
coastline. Since significant waveheight is typically depressed from 
0.5 to 1.0 meters in the oil spill region, the measurement of signi- 
ficant waveheight is closely coupled with the analysis of oil 
spills. The Coast Guard is currently developing a three-dimensional 
model for oil spill distribution, and historical significant 
waveheight data are used as an input to that model. 
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Naval architects, sea engineers and ship builders also 
need historical waveheight data. Since ships and ocean structures 
are designed to withstand the structural stresses that they might 
reasonably be expected to encounter, a history of waveheight by geo- 
graphical area and season could be quite helpful. A OEOS-3 atlas of 
SWH is believed to be the most accurate source of historical 
waveheight data currently available and is, therefore, the best 
source for design information. Although subtle differences between 
the GEOS-3 significant waveheight atlas and other SWH data bases 
might not seem to be of sufficient magnitude to affect the design of 
ships and offshore structures, the opposite may be true. Dr. Herb 
Austin of the Virginia Institute of Marine Science has stated that 
the OEOS-3 atlas would be ^extremely useful” in this context (Aus- 
tin, 1979). 

As previously mentioned, some areas of the world’s oceans 
are so remote that significant waveheight data from these regions is 
virtually nonexistent. No ocean buoys are located in these areas 
and ships rarely traverse through them. The U. S. Coast Guard has 
expressed interest in analysing significant waveheight contour data 
from these remote areas, with particular interest in how meteorolog- 
ical features move through these rarely traveled regions. 

Oceanographers, meteorologists, and climatologists fre- 
quently plan their research, field trips, etc., based upon the sta- 
tistical probability of locating certain climatic conditions. The 
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GEOS-3 atlas could be used by scientists to plan research when pred- 
iction of significant waveheight is important. 

Finally, the GEOS-3 near-real-time significant waveheight 
estimate has already gained wide acceptance in the scientific com- 
munity as an accurate and timely data set to which other significant 
waveheight measurements can be compared. The GEOS-3 historical sig- 
nificant waveheight contour atlas presented in this report could 
similarly be used as a reference point for comparison of other 
results. 

1 . 7 Scope of the Investigation 

The first part of this study attempts to document the 
accuracy of the GEOS-3 altimeter derived significant waveheight 
estimates. The estimates, which were obtained from a statii&tically 
representative sampling of GEOS-3 passes, were compared with 
independent measurements of significant waveheight. These indepen- 
dent measurements were obtained from the NOAA Data Buoy Office (Had- 
sell, 1974) which f.i&intains e number of ocean stations (see Figure 
1.4) that routinely measure significant waveheight and other sea 
state and atmospheric purametera at 3-hour intervals. The accuracy 
of the buoy measurements of sigr.ificant waveheight has been esta- 
blished to be 55 cm (Steele and Johnson, 1977). 

For the purposes of this pai' i of the investigation, it was 
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necessary to identify 


1. all GEOS-3 passes whose ground tracks came arbitrarily close to 
one of the NOAA buoy locations, and 

2. those passes identified in (1.) where the altimeter was track- 
ing in the intensive mode. 

It was decided that in order to satisfy both of the above conditions 
the GEOS-3 data must have been taken within 90 minutes of a buoy 
measurement and must have passed within one equatorial degree (111 
kilometers) of the buoy making the measurement (see Section 3«1)* 
The estimates of significant waveheight from these GEOS-3 passes and 
the corresponding buoy measurement of significant waveheight form 
the data set for the first part of the investigation. 

Finally, the GEOS-3 significant waveheight estimates for 
the entire mission are presented in the form of a global atlas. 
This atlas, patterned after the significant waveheight results 
illustrated in the U. S. Navy Marine Climatic Atlas (U. S. Navy, 
1974), presents global significant waveheight contour maps for both 
low and high sea state. 
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CHAPTER 2 

SIGNIFrCANT WAVEHEIGMT FROM SATELLITE ALTIMETERS 
2. 1 Wavcforw Geometr y for Negl igtble SW| 

The geometry of a square |>ulse emitted from an altimeter 
antenna which lm^>ingea upon an idealized flat sea surface is illus- 
trated in Figure 2.1. This figure depicts the distances from the 
spacecraft to the subsatellite point (the point on the surface of 
the earth lying on the line between the spacecraft and the center of 
mass of the earth) and to a general point P near the subsatellite 
point . 


In Figure 2.2, a pulse of duration t is observed leaving 
the spacecraft at time t = 0 and returning from the sea surface at 
time t = tg , after traversing a distance H each way. At the 
surface, it is reflected (at time t = t^ ) back to the satellite 
where it is received at the antenna. Assuming that any vertical 
motion of the spacecraft has negligible effect upen the signal tran- 
sit time, 

tp » ^ (2.1) 

c o 

where c is the speed of light. At time t = t^ = t^ 4 - iS , where 
6 Is some time increment less than T , the square pulse is 
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FIGURE 2.2. tOEALIZEO MEAN RETURN PULSE SHAPE 



observed impinging upon the see surfaoe. The observed area 
iUuminated by the pulse expands ulroularly until the trailing edge 
or the waveform Is received at the satellite. If the distance from 
the spacecraft to point P on the edge of the llluminited circle 
(see Figure 2*1) is H>h, then the radius r of the illuminated 
area is related to H by 

r^ s (H 4 h)^ - : 2Hh (2*2) 
since h << H and therefore h^ « 2Hh . 


Tl\e time t = t^ is given by 


. . 2H 

t.j s t« '= -f 

3 2 o 


( 2 * 3 ) 


and corresponds to the two-way travel time between the satellite and 
point P on the surfaoe. TJ^erefore, 


1 ,2H .t* H ♦ h 

2 o c 


( 2 * 4 ) 


or 


M H h 

c 2 o o 


which yields 



Combining Equations (2.2) and (2*6) yields 


( 2 . 5 ) 


( 2 * 6 ) 


an 


r^ s Ho 6 


( 2 * 7 ) 



and the area of the illuminated region is given by 

A a ^r^ a nHo^ (2.8) 

Note that the observed area of the illuminated region, or 
equivalently the power received by the satellite, increases linearly 
with time. 

At t tj| , the trailing edge of the pulse is received at 
the satellite. At time t a t^ where t^ > tjj , the observed 

illuminated area becomes an annulus with inner radius rj end outer 
radius r^ given by 

Tj a Hc(t^ - tj^) (2.9) 

rl a Hcit^ - t^) (2.10) 

and width W and area A given by 

W a rg - rj = (Ho(t5-t2)]'^^ - [Hoit^-ti,)] (2.11) 

A a nr^ - Tir^ a ■nHcT (2.12) 

From Equation (2.12), the area of the annulus remains constant. 
Therefore, the power received at the satellite remains constant 
until antenna beamwidth effects cause the power received to decay. 
This is due to the increasing size of the annulus and the limited 
antenne beamwidth (which for GEOS-3 was 2.3^). 
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The character of a square pulse impinging upon a sea sur- 
face with negligible significant waveheight as seen from the satel- 
lite can be summarized as follows: 

1. no power is received until the leading edge of the pulse 

strikes the sea surface and is reflected back to the satellite; 

Z. after the leading edge of the pulse is reoeiveu at the space- 
craft and before the trailing edge is, the power received 
Increases linearly with time; 

3* after the trailing edge of the pulse is received at the space- 
craft, the power received remains at a constant plateau value; 

and 

after the antenna beamwidth effects become non-negligible, the 
power received begins to decay. 

These four stages are depicted in Figure 2.2, which represents the 
idealized mean return pulse shape. However, due to the scattering 
properties of the surface, the Instantaneous power received fluctu”' 
ates, making it necessary to average a large number of pulses in 
order to determine the mean pulse shape. 

The GEOS-3 satellite received the return pulses in 16 

waveform sampling gates (see Figure 1.3). These 16 values, called 
Instantaneous Return Samples, or IRS’s, were collected 100 times per 
second by the spacecraft and averaged onboard in an attempt to 
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construct an accurate representation of the mean ipulse shape, which 
was characterized by the four stages described in the previous para- 
graph. The 16 values of averaged IRS's were called Average Return 
Samples, or ARS's, and were computed using an RC filter with a 2- 
second time constant (teitao and Purdy, 1975). It will be shown in 
Section 2.2 that the significant waveheight can be determined by 
observing the departure of the leading edge slope of the ARS's from 
the leading edge slope that the ARS's would have in the ideal calm 
sea case. 

2.2 Waveform Geometry for Non - Negligible SWH 

In the case where a square pulse impinges upon an ocean 
surface with non-negligible significant waveheight, the shape of the 
mean return pulse will be altered. The geometry of a square pulse 
impinging upon a sea surface with non-negligible significant 
waveheight is illustrated in Figure 2.3. Note that the crests of 
the waves are illuminated prior to the time at which the calm sea 
would have been illuminated. Similarly, the troughs of the waves 
are not illuminated until after the time at which the calm sea would 
have been illuminated. The net result of these effects is that the 
mean power received for non-negligible sea state does not reach its 
full plateau value until after the time at which the mean power 
received from a calm sea would have reached its plateau value. This 
effect causes the slope of the leading edge of the ARS’s to dimin- 
ish. 
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The leading edge of an idealized mean pulse shape for 
negligible sea state and for several non-negliglble sea states is 
characterized in Figure 2 ,^, From this figure it can be seen that 
the slope of the mean return pulse is related to significant 
waveheight. If the mean pulse shape for negligible sea state is 
knovm precisely, then the significant waveheight can be determined 
by analyzing the departure of the mean return pulse shape for the 
non-negligible sea state from the mean return pulse shape for calm 
sea. 

2.3 Choice of Mode l 

It has been shown (Brown, 1977) that the mean return 
waveform can be conceived of as a convolution of 

1. the system point-target impulse response, 

2. the non-coherent surface impulse response, 

3. the Ocean surface height probability density function, and 

4. the tracking loop Jitter. 

The first component is the composite of the transmitted pulse and 
the transmitter and receiver bandwidth effects. Its distribution is 
a complicated function which can be modeled as a Gaussian distribu- 
tion. The second component is the calm sea response, which can be 
modeled as a step function. The third component represents the 
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rough sea distribution, which is modeled in this investigation as a 
Gaussian distribution. The fourth component is assumed to be unaf- 
fected by sea state in this investigation. (See Section 2.6 for a 
discussion of the errors associated with assuming these distribu- 
tions.) 


Assuming that pointing angle errors have negligible effect 
upon the leading edge of the waveform (see Section 2.6), Brown and 
Miller (197^) have shown that a good approximation for the return 
power as a function of time is the Gaussian function 

y(t) = aP(-”) ♦ d (2.13) 

where 

a = return waveform amplitude 
b s time origin 

c = return waveform standard deviation 
d = return waveform baseline amplitude 
and where p(^^) is the Probability Integral 
z 

P(z) = / Z(q)dq (2.14) 

-KS 

and Z(q) is the Gaussian function 

Z(q) = e"*’ n (2.15) 

(27T)^'^ 

P(z) is determined (from Abramowitz and Segun, 1968) by 


33 



where 


PU) s 1 


- Z(z)(b^t 


“>2* 


bjf 


l>„t 


bjt*) 


♦ E(l) (2. 16) 


- JL 

(2.17) 

Up 

c(z)t < 7.5x10"® 

(2.18) 

= 0.2316419 

(2.19) 

y s 0.319381530 

(2.20) 

2 = -0.356563782 

(2.21) 

= 1 .781477937 

(2.22) 

= -1.821255978 

(2.23) 

>5 = 1.330274429 

(2.24) 


In this model it Is necessary to estimate the four parame- 
ters a, b, c and d from which the significant waveheight can be 
determined. The technique used to estimate the four -parameter func- 
tion y<t) is the method of least squares. 


2**l Derivation 


Equation (2.13) can be expanded in a first-order Taylor 

Series approximation about a point y = y(a ,b^,c ,d ) (see Hayne, 

o o o o o 

1977) to yield 


^i 


y ♦ y (a-a ) y^(b-b ) y (c-c ) 
0 a o b 0 0 0 


♦ y*(d-d ) ♦ . . • 
'd o 


(2.25) 
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where 


■i > P(-r-2) 

y=y„ 


y s 

8b 


•ft ^”^ft 

r 

o o 


8 c 


a t-b t-b 

^ ( 2) z 

®o % ®o 


y s 

8d 


and where are the observed values of y , y^ are the computed 
values of y , and higher order terms have been neglected, thereby 
assuming that y^ is sufficiently close to y^ to permit ^jonver- 
gence while neglecting the complications introduced by including 
higher order terms in the Taylor Series expansion, The traditic lal 
least squares estimate of y(t) is obtained by minimizing the sum 
of the errors squared over the 16 gates 


E * 2 (y. - y,)‘ 

i=1 ^ ^ 


(2.30) 


with respect to (a-a ), (b-b ), (c-c ) and (d-d ) . When this is 

o o o o 

done, the following four equations are generated: 


= 0 = 22 (y,-y,)-; 


8(y.-y.) 16 _ 




= 2 2 


(2.3t) 
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^ 

3(b-0 


^6 _ 

1 s 1 O 1 s 1 


(2.32) 


3E 


16 


^(yryt'f - 

I 1 _ ^ 


5f^ = 0 . g^z^(yr>i>-3T5--'i;r ' ‘^- 


33) 


3E 


16 __ 9(y,-y.) 16 

s 0 s 2 2 (y -y ) y - s 2 2 (y.-y. )y 

is1 ' ' 1 1 


(2.3*0 


Substituting Equations (2.26) through (2.29) into Equations (2.31) 
through (2.3*1) yields 


16 


16 


16 




16 


16 


16 


♦ 2 y-yQ^c-c ) + 2 ygyd^^'^^o^ ” ^ ^i^a 

i_1 a c o i-1 a a ° i a 


16 


16 


16 


° = 2 y.yj, ♦ 2 yayu^a-a ) ♦ 2 y. (b-b ) 


16 


16 


16 


♦ 2 ywy«<c-c ) ♦ 2 - 2 y.y. (2.36) 

bo o b d o i b 


16 


16 


16 


° = 2 y„y. + 2 y y (a-a ) •»■ 2 ywy.(»>-b ) 

0 c a c o bo o 


16 2 16 

+ 2 yc<c“Co> + 2 ycyd^'‘"‘*o^ " 2 yiyc 
i=l ° ° 1=1°** ° 1=1 ^ ° 


( 2 . 37 ) 


16 


16 


16 


° = 2 y„y^ •»• 2 y«yH(a-a«) ♦ 2 ywyn^^^'^o^ 

j_^ad o ^_^bd o 
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(2.38) 


16 


16 


16 


is1 1*1 1*1 


Equrttions i?,i3b) through (2.38) cnn be written in matrix form as 


16 2 16 16 16 

2, Vc 

l*«l i»l l«l 1*"1 


16 


16 


16 2 

,’\Vb ,’\Vc 

1-1 1-1 1-1 i-l 


16 


16 


1 6 2 ^ 6 
i^l ** C lm\ *’ 1-1 im\ “ 


16 


16 


16 


16 


/, V.1 5'c’'d 

i-l i-l i-l i-l 


a- a 


b-1) 


C“C 


d-d 


16 

l-l 

16 

i-l ^ ^ ° 

16 

yc^yryo> 

yd<yryo> 

1-1 


( 2 . 39 ) 


For thin problem* the technique of least squares is 

applie^d to a tiHmoated Taylor series (Equation 2.25), which is 

linear in the corrections to the guessed values of a ♦ b , c at>d 

o 0 0 

d^ and, therefore, separable. Aooordinftlyt iteration is required 
because of the 1 inear i nation. 


It should be noted that in order to solve for a, b, o and 
d, a 4 X 9 symmetric matrix must be inverted, Tliis matrix has been 
examined and found to be well conditioned. In particular, the accu- 
racy of the matrix inversion has been examined for a wide range of 
waveheights and was determined to be satisfactory. 

Tl^e convolutional model for the mean return waveform (see 
Section 2.3) assumes that the standard deviation, c (sometimes 
referred to ns the risetime coefficient which should not be confused 
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with th<* spt’od of light)# is a composite of four elements which can 
be grouped in the following way: 

1. the ocean surface height probability density function, and 

2. the composite of a) the transmitter and receiver effects# b) 
the noncoherent impulse response and o) the tracking loop 
J i tter , 


Tile first of these two groups can be referred to as the ’’rough sea” 
contribution to the return power and the second as the ’’calm sea” 
contribution to the return power* Since both distributions are 
assumed to be Gaussian and since the convolution of two Gaussian 
distributions is Itself a Gaussian distribution, 

? of ♦ (2.110) 

s 0 

where o is the calm seo standard deviation, and o is the rough 

0 s 

sea standard deviation, both expressed in nanoseconds. 

According to Neumann and Pierson (1966) 


The significant waveheight is defined as the average of the 
heights of the one-third highest waves In a long seQuence of 
Waves observed at a point. It is more or less equal to four 
times the square root of the variance of the wave record. 
Tests of wind-sea records, by first averaging the heights of 
the one-third highest waves and then. estimating m (the wave 
record variance) and computing 4 [m 3 , yield th^ jjame value 
within perhaps 5 per cent. 

Using Equation (2.i|0) and converting to units of meters by multiply- 
ing by the two-way speed of light yields: 
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SWH s 4 0 * O.b tc’^- 

Thtj iinplomv»'nt«tlon of the estimation algvn'ithm proceeds as: 


1. 

Provide 

initial 

guesses for a , b , 

0 o 0 

and d 

0 

2, 

Estimate 

a. 

b, 

c and 

d using Kquation (2 

.39) 

3* 

Compute 

E 

usl 

ng Equation (2.30) 


4. 

If E has c 

on verged , 

go to Step #6 



Hepl ace 

a 

o 

, b 

0 • ^0 

and d with the new 
0 

estimates of 


c and d and return to Step 12 
b. Compute SWM using Kquation (2,41). 


2, 5 Converj^ence C^onsiderat i tu 

Since Equation (2,2^) negleoted higher order terms, the 
estimation algorithm must be iterated, Tlie oonvergence criteria for 
the estimation algoritiim Uvas that the relative error in Equation 
(2. .10) for two oonseoutive iteratives be less than O.lt, Signifi- 
cant waveheight has been computed in this way on thousands of passes 
of (iEOS-3 data, and it has been found that the algorithm nearly 
always convergos within 2 or 3 iterations. In addition, it has been 
determined that the final converged estimate of significant 
waveheight is not particularly sensitive to the initial guesses of 
**o ' ^o ’ % ^o ' practice, the following initial guesses 



are for all passes: 


a s 84. *500 (2.42) 

o 

s - 0,902 (2.43) 

2 8.500 (2.44) 

d s 5.800 (2.45) 

o 


For each significant wavohelght estimate after the first frame of 

data, the converged values of a, b, c and d for the previous frame 

of data are used as the initial guesses for a , b , c and d 

0 0 0 o 

As was discussed in Section 2.2, it is necessary to have 
an accurate estimate of the calm sea standard deviation, , in 
order to calculate significant waveheight using Equation (2.41), 
Early in the GEOS-3 mission, many arcs of the satellite which passed 
over areas where ship measurements indicated the presence of calm 
seas were analyzed in order to determine that value of 0^ which 
would yield an estimate of SWH = 0 for those passes. The value 
arrived at was 

a =7.49 nsec (2,46) 

o 

This value, which accounts for the effects of tracking loop jitter 
in an average sense, was examined by McMillan and Roy (1977) and 
found to produce bettor agreement with ground truth measurements 
than did any other tested value. Gower (1979) estimated that 



lomoving the olTects of tr^icking loop jitter would yield 


0 s 6 . 05 + 0.2b nsec (2.H7) 

c ■ “• 


Rocause rough sea scattering behaves as if it were a col- 
lection of incoherent, discrete scntterers and the .statistical pro- 
perties of the scattering ar'o assumed to be Rayleigh (see, for exam- 
ple, Walsh, the mean return represented by the AR$*s contains 
noise. Although individual waveforms exhibit these Rayleigh fluc- 
tuations, the data points are all scattered about the mean return 
for some averaging interval (Hammond, et al, 1977). Nevertheless, 
the waveforms from two adjacent data records of ARS's can differ 
substantially, even though the altimeter is receiving data from 
ocean areas only a few Kilometers apart, lliis, combined with the 
numerical errors associated with the estimation process, produce an 
estimate of c which sometimes causes the term under the radical in 
Equation (2.41) to become negative if the true sea state is very 
calm , 


Figure 2.5 illustrates the algebraic relationship between 
significant wavehelght and the estimated value of o . It should be 
noted that for moderate to large values of significant wavehelght, 
small errors in the estimate of c do not cause large errors In the 
calculated value of significant wavehelght. However, the estimate 
of aignificant waveheight is very sensitive to even small errors In 
c for calms seas, i.e., for significant waveheight less than about 1 
me t e r , 


4.1 
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When the scattering properties of the .surface, the noise 
characteristics of the altimeter, the estimation errors and the 
algebraic sensitivity of the estimate to c are combined, it is 
obvious that smoothing the estimate should provide more confidence 
in its accuracy, especially when the sea state is calm. The longest 
segment over which the sea state can be assumed to be correlated has 
been empirically determined (Apel, 1975) to be 150 kilometers (or 21 
seconds for GEOS-3). Therefore, a sliding 21-second rectangular 
filter was employod by the significant waveheight estimation 
software. 


Either the estimate of c or the calculated value of SWH 
can be smoothed and the results were shown by McMillan and Roy 
(1977) to be identical to within the numerical precision of the com- 
puter. For computational ease, the estimate of c is smoothed in 
this investigation. However, even when the estimate is smoothed, 
the term under the radical in Equation (2.1<1) .still occasionally 
become.s negative due to the effects of noise in the measured 
waveforms and the algebraic sensitivity demonstrated previously in 
this section. Such estimates had no physical meaning and a value of 
SWH = 0 was assumed. 

The first few weeks after the launch of GEOS-3 were desig- 
nated the calibration phase of the mission and were designed to 


43 



eliminate known errors and inconsistenf' a in the preprocessing of 
the altimeter data. During this phase, it was determined that two 
important corrections were required for the ARS's. 

First, due to limitations of current circuitry and to tim- 
ing difficulties, the 16 sample and hold gates should not have been 
assumed to be equally spaced in time. Accordingly, General Electric 
supplied WFC with the timing corrections necessary to properly time 
tag the ARS's (see Table 2.1) based upon correlation analysis of the 
ARS's. 


Second, it became evident that amplitude biases needed to 
be determined for the ARS's. In subsequent weeks, several sets of 
ARS amplitude biases were determined by G. S. Hayne and by E. J. 
Walsh at WFC (Walsh, March 1979). McMillan and Roy (1977) deter- 
mined that the most consistent agreement between the estimated value 
of significant waveheight and direct measurements of significant 
waveheight made by buoys and ships was produced when the Walsh 
amplitude biases were employed (see Table 2.2). Since the ship and 
buoy measurements of significant waveheight constituted a statisti- 
cally representative sampling population, it was assumed that "tun- 
ing" the estimation algorithm to that particular set of ground truth 
data would not cause an estimation bias. The validity of this 
assumption was later demonstrated when several independent estima- 
tion algorithms were compared (Fedor, et al, 1979) and were found to 
produce essentially equivalent estimates. 
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ARS GATE 

TIME (nt,^ 

ARS GATE 

TIME (ni) 

1 

-52.19 

9 

- 6.88 

2 

-46.00 

10 

^ 0.00 

3 

-43.63 

11 

^ 6.50 

4 

-37.50 

12 

f12.09 

5 

-31.81 

13 

4^15,19 

6 

-24.88 

14 

425,69 

7 

-17.12 

15 

431.69 

B 

-12,31 

16 

438.38 


TABLE 2.1, ARS RELATIVE TIMES 


ARS GATE 

BIAS imv) 

ARS GATE 

BIAS (mv) 

1 

4 2.3 

9 

4 1.3 

2 

- 2.7 

10 

- 2.0 

3 

4 0.8 

11 

4 3.6 

4 

- 1.8 

12 

4 1.3 

5 

4 2.5 

13 

4 0.9 

6 

- 0.1 

14 

- 0.5 

7 

-0.8 

15 

- 0.3 

8 

- 1.2 

16 

- 4.0 


TABLE 2.2. ARS AMPLITUDE BIASES 



The algorithm derived in this chapter, together with the 
smoothing technique and ARS timing and amplitude bias corrections 
detailed above, was used at the Wallops Island facility in prepro- 
cessing the GEOS-3 altimeter data (McMillan, i975). The same tech- 
niques were used by the Goddard Space Flight Center in the near- 
real-time data network (McMillan, 1978) described in Chapter 1, 

2,6 Error Sources 

The error sources affecting the estimation of significant 
waveheight can be divided into three general categories: 

1. modeling errors associated with the waveform shape 

2. measurement errors associated with the return power 

3. mathematical errors associated with the estimation of the slope 
of the ARS's 

Each of these categories will be examined separately. 

The first category of error sources includes the errors 
associated with the modeling of the waveform shape. As presented in 
Section 2.3, the mean return waveform is modeled as a convolution of 

1. the system point-target response, 

2. the noncoherent surface impulse response. 


46 



J. the ocean surface height probability density function, and 
the tracking loop Jitter. 

The first of these terms, which represents a composite of the 
transmitted pulse and the transmitter and receiver bandwidth 
effects, resembles a Gaussian distribution and is therefore modeled 
as such, introducing a model error. The second term is the calm sea 
impulse response, which resembles a step function and is so modeled, 
introducing another error. The third term is the radar observed 
distribution of the ocean surface. The distribution is assumed to 
be Gaussian in nature. If the distribution is not Gaussian, i.e., 
if the probability density function must be characterized by higher 
order terms (skewness and kurtosis), then the estimation process 
becomes more complicated. The skewness has been accounted for by 
Walsh 0979). Nevertheless, Fedor, et al ( 1979) found that the 
Walsh algorithm and the Wallops algorithm produced nearly identical 
estimates. For a further discussion of the surface elevation proba- 
bility distribution, see Huang and Long (1981). 

Additionally, the radar observed ocean distribution is not 
the true geometrical distribution. This difference might be 
accounted for through the use of a correction called the electromag- 
netic bias (EM bias), which is currently being investigated (Jack- 
son, 1979) for use in reducing the altimeter data of future space- 
craft. The effect of the EM bias is that (McMillan, et al , 1980), 



th« relative radar croas section tends to increase below mean 
sea level and decrease above mean sea level in the presence of 
waves. Its effect is to shift the centroid of ^he radar return 
away from mean sea level toward the wave troughs so that the 
altimeter tracks long. Recent experimental data from the Sur- 
face Contour Radar at 36 GHz and the NRL 10 GHz adaptive radar 
altimeter itidicate that the EM bias Is in the range of 0 to 3 
percent of the SWH. 

Each of the first three components in the convolutional 
model of the return waveform is assumed to have a defined functional 
form which can easily be incorporated into an overall model of the 
waveform. Actually, each of these modeled functional forms intro- 
duces an error into the estimation process, but a more accurate 
model for the estimation of significant waveheight from GEOS-3 data 
remains to be established. 

As discussed in Section 2.3« the tracking loop jitter is 

assumed to be independent of sea state and is not accounted for 

directly in the GEOS-3 significant waveheight algorithm. It is 

accounted for, in an average sense, in the determination of the calm 

sea pulse width, <T (see Section 2.5). Any error introduced in 

0 

the estimate by not properly accounting for tracking loop jitter 
would tend to be more significant for high sea states than for low 
sea states (Hayne, 1976). Walsh (1979) and others have eliminated 
the effects of tracking loop jitter in their SWH estimation algo- 
rithms. Fedor, et al (1979) examined both actual and simulated 
GEOS-3 data and determined that this additional numerical procedure 
did not appreciably alter the estimate (see Section 2.8). 
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Additionally, tht* Gaussian distribution used to model the 
return power is presented as a function of four parameters, each of 
which must be determined In the least squares estl nation process*, 
Other Gaussian functions with a different number of unknowns could 
be used bo model the return waveform, bub one would expect that no 
significant variation in the estimated SWH would result by altering 
the algorithm In this manner. 

The second general category of error sources includes the 
errors associated with the measurement of the return power. In Sec- 
tion 2,5, the ARS amplitude and timing biases were presented. 
McMillan and Roy (1977) examined several sets of ARS amplitude and 
timing biases and concluded that results varied by as much as 15< 
when using different bias sets. Both of these bias corrections to 
the measured return waveform are somewhat arbitrary, especially the 
amplitude biases. Nevertheless, these corrections are the most 
accurate biases currently known. 

Certain other instrument-related error sources can be 
enumerated under the general category of errors associated with the 
measurement of the return waveform. They include the effects of 
pointing angle, pulse width, gate position, AGC fluctuation, power 
and timing variations and gate saturation. For GEOS-3, the 
estimated standard deviation, c , was assumed to be insensitive to 
variations in pointing angle. The validity of this assumption has 
been proven by G. S. Hayne at WFC using simulated waveform data with 



varying pointing angles. The pulse width is a limiting factor in 
the resolution of the estimate, with a smaller pulse width yielding 
higher resolution. The position of the sample and hold gates (see 
Figure 1.3) was not optimum for the estimation of significant 
waveheight, and therefore is a potential error source. A more 
optimum gate configuration would include more sample and hold gates, 
especially in the ramp portion of the waveform, which is the portion 
most sensitive to sea state. AGC fluctuation, power and timing 
variations and gate saturation could and did occur at certain times 
during the mission, but these effects were assumed to cause only 
noise in the SWH estimation algorithm (Stanley, 1980). 

The third and final category of error sources includes the 
mathematical errors associated with the significant waveheight esti- 
mation process itself. These include all numerical estimation 
errors such as matrix inversion errors, roundoff errors and trunca- 
tion errors. All of these errors are assumed negligible. The vali- 
dity of this assumption is supported by the facts that the Wallops 
algorithm determined equivalent SWH estimates on at least three dif- 
ferent computers, and the Wallops algorithm compared favorably with 
the other estimation algorithms examined by Fedor, et al (1979). 
Additionally, the accuracy of the matrix Inversion has been substan- 
tiated over a variety of ciea state conditions. Finally, it should 
be pointed out that G. S. Hayne at WFC is currently Investigating 
the contributions of the various error sources in the significant 
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Tlu’ prublam of cnUnilniing oa<nvn significant wavcheight 
hastHl upon the altimeter waveform measurements of the yE0S~3 space- 
craft has been addressed by a number of other GEOS-3 principal 
Invest i gator s (Kedor and nat'rie.k, 1978; Gower* 1979; Hayne, 1977; 
Rufenaeii and Alpers, 1979; and Vialsh, duly 1979). 'Hre various tech- 
niques used by the GKOSOj principal investigators were reviewed dur- 
ing the design of the Seasat altimeter t As a consequence* the sig- 
nificant wavelrelght was calculated onboard the Seasat spacsccraft 
using the altimeter return waveforms in much the same fashion as has 
been presented in this study. Although the purpose of this investi- 
gation is not to uompare the accuracy of various algorithms but to 
establish the accuracy of the WEC and real-time significant 
waveheight algoi'ithm* n description of tiiese other algorithms is 
included here for the sake of completeness. A more detailed compar- 
ison of the algoritiwis is given by Fedor, efc al , (1979). 

Although tl»e different algorithms solve for different sets 
of parameters, use different weighting and best fit criteria and 
even dvffer as to whether IRS or ARS data are used as input, they 
all fit a model function to the detected waveform. Tlius, the 
(Ufferences in tlie accuracy of the algorithms are essentially due to 



the curve fitting techniques employed. Accordingly, the comparison 
study (Fedor, et al, 1979) found that "individual differences 
between the algorithms were small when compared to the general good 
agreement among them." 

As previously noted, several of the algorithms employ the 
ARS's as the input to the curve fit technique, while the other algo- 
rithms use the IRS's. For 3.2 seconds of data, a fit through the 
ARS data requires the determination of a curve passing through only 
16 points, whereas tha use of IRS data requires the determination of 
a function which best represents 5120 (320 x 16) points. Not only 
does the use of IRS data slow down the estimation process due to the 
requirements of reading and storing so many variables and burdening 
the curve fit software with such a cumbersome number of data points, 
it precludes the use of telemetry mode data (low data rate data) 
which does contain ARS's but does not contain IRS's. 

The GEOS-3 altimeter mode could be 

1. Global 

2. Intensive 

a. none of the IRS's reported 

b. half of the IRS's reported 

c. all of the IRS's reported 

If the altimeter were operating in the global mode, no waveform 
Information was reported so that significant waveheight could not be 
estimated. If the altimeter was operating in intensive mode, ARS's 


WiM e iilwuys IRS's may or may not have been repotted. 
Thus, the algorltlims that processed IRS data were not only limited 
to intensive itigh data rate but were al.so limited to those intensive 
high data rate modes wltere all of the IRS's were reported. 

It sitould also be noted that the different algorithms 
treat the effects of tracking loop )Utt*r in different ways. 
Although tile WFC sigriificant waveheight algorithm (as well as 
several of the other estimation algorithms) treat tracking loop 
jitter in an average sense by including its effects in the calm sea 
risetime (see Kquation 2J\6) and assuming that it is unaffected by 
the magnitude of significant waveheight, several of the algorithms, 
i.e., Walsh (July, 1979), have attempted to acc' \it for the contri- 
bution of tracking loop Jitter in a more rigorous way. These tech- 
niques involve a realignment of the IRS’s based upon the residuals 
between the actual altitude measurements and smoothed altitude meas- 
urements. Since they use IRS data, the 'FM Mode 2 (high data rate) 
must be used. Nevertheless, the close agreement between these 
methods of estimation and the methods which accounted for Jitter in 
an average sense seams to indicate that Jitter need not be deter- 
mined rigorously. 

2.8 Coin^arison of Other S^ Esti mation Techniques 

Of the six significant waveheight estimation algorithms 
presented by Fedor, et al (1979), three of the algorithms solve for 



four or more parameters by fitting a curve to the leading edge of 
the return waveform. Four of the algorithms use a least squares 
technique to achieve convergence of the curve fit, and four of the 
six employ the IRS*s instead of the ARS's, although one of those 
models uses an abbreviated set of IRS's. The major characteristics 
of each of the algorithms is summarized in the following: 

Wallops Algorithm . The WFC significant waveheight algorithm 
presented previously in this chapter employs a Jesot squares 
fit of a four-parameter function to the ARS data. Tracking 
loop Jitter is accounted for in an average sense. This algo- 
rithm Was used for the computatioh of SWH for all archived 
GE03-3 data and for the real-time estimation of SWH. 

Walsh Algorithm (Walsh, July 1979)* The Walsh significant 
waveheight algorithm performs a five-parameter least squares 
fit to the IRS data. Tracking loop Jitter is accounted for by 
a time realignment of the IRS's based upon the residuals 
between the actual altitude measurements and smoothed altitude 
measurements . 

Gower Algorithm (Gower, 1979). Hie Gower significant 
Waveheight algorithm performs a modified maximum likelihood 
estimation of a four-parameter fit to the IRS data. Time 
rualignment of the IRS's is used to account for tracking loop 
Jitter. 



Rufenach and Alpers Al gorithm (Rufenbach and Alpers, 1979). 
The Rufenach and Alpers algorithm performs a least squares 
estimation of a two-parameter fit to IRS 8 through IBS 12, 
which Is the rise-time portion of the return waveform. 
Realignment is performed to account for tracking loop jitter. 

Fedor Algorithm (Fedor and Barrick, 1978). The Fedor signifi- 
cant waveheight algorithm fits a three-parameter function to 
the first differences of the IRS's after a time-realignment 
correction for tracking loop jitter. An iterative procedure 
then computes a correction to account for the difference 
between first differences and true point derivatives. The 
best-fit function is determined by the method of least squares. 

Godbey Algorithm . The Godbey significant waveheight algorithm 
is based upon a monotonlcally increasing function of the ARS's 
and thus accounts for tracking loop jitter in only the average 
sense. The algorithm estimates only one parameter and does not 
require iteration. This computational simplicity is a real 
advantage, especially over the algorithms Which process IRS 
data. However, some of the contributors to the Fedor, et al 
(1979) investigation (Hayne, I960) believed that the Godbey 
algorithm simplified the estimate too much and that the actual 
shape of the waveform should be estimated. 


Fedor, et al (1979) concluded that all of the algorithms 



agreed well with each other and with a small amount of ground truth 
data which was available for their evaluation. The study esta- 
blished the following standard deviations for all of the algorithms; 

a < 0.75 m. for 0.0 < SWH < ^.0 m (2.48) 

a < 0.50 m, for=4.0 < SWH < 8.0 m (2.49) 

From the results presented in this chapter and later in 
Chapter 4, it can be deduced that since the significant waveheight 
estimated by the Wallops algorithm agrees well with NOAA buoy data 
and since the six algorithms presented in the above study all agree 
with each other, that all of the algorithms are providing satisfac- 
tory estimation accuracy. Therefore, simplicity of operation 
becomes an overriding concern. This would seem to indicate that the 
algorithms which process ARS data are not only preferable because 
they allow processing of low data rate data, but also because they 
achieve comparable accuracy without the computer time and storage 
requirements necessary for processing IRS data. 
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CHAPTER 3 


THE DATA SET 


3. 1 GEOS Data Set 

For the purpose of estimating the SWH measurement accu- 
racy, it was necessary to select those passes of GEOS-3 altimeter 
data from the entire set of GEOS-3 data segments (see Figure 3.1) 
which could be compared to the independent buoy measurements of sig- 
nificant waveheight. The GEOS-3 ground tracks, which were computed 
at GSFC to a radial accuracy of 1 to 1.5 meters (Lerch, et al, 
1978), must have passed near one of the buoy locations (see Section 
3,2) at about the time when the buoy waveheight measurements were 
made . 

As was discussed in Section 1.^, the GEOS-3 altimeter did 
not operate continuously but was subjected to the power constraints 
of the spacecraft. Accordingly, the first criterion for matching a 
GEOS-3 data pass with a buoy observation was that the altimeter be 
operational at the time that the ground track passed the buoy. In 
addition, the altimeter must have been operating in a mode which 
allowed significant waveheight to be estimated. 

Since the magnitude of wind generated waves in the open 
ocean varies slowly with respect to time and distance, it can be 


57 




58 



nasumed to be unchanged over areas of moderate size surrounding the 
buoys. Therefore, the radius of the search area about the buoy 
location Is somewhat arbitrary. A search area corresponding to one 
equatorial degree (about 111 kilometers) in radius was assumed to be 
sufficiently small such that the significant waveheight was 
unchanged throughout the entire area under reasonably normal condi- 
tions. The validity of this assumption, which has been verified by 
Apel (1975), will be considered in Section 3.2. Additionally, a 
search area with a radius of one degree is large enough to identify 
a statistically representative (but not cumbersome) sample popula- 
tion of altimeter data passes. 

The mere passing of the satellite ground track through any 
of the buoy search areas did not qualify that pass as a member of 
the comparison data set. Since the buoys were often irregular in 
measuring significant waveheight and since all of the buoy measure- 
ments during the GEOS-3 mission were not readily available (see Sec- 
tion 3.2), it was necessary to identify those altimeter !ata seg- 
ments passing through the search area at a time when the buoy was 
measuring significant waveheight rnd at a time when the buoy meas- 
urement was available for reduction. 

If the altimeter data set was restricted to data segments 
whose ground tracks passed through the search area at the exact time 
when the buoy measurement was being made, there would be only a 
small number of passes to be reduced. Fortunately, waveheight is 
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normally a slowly varying parameter in the absence of storm condi- 
tions. Therefore, any altimeter data segment which passed through 
the search area and which occurred within some reasonable period of 
the buoy measurement could be accepted in the altimeter comparison 
data set. This period was chosen to be 90 minutes, since the buoy 
SWH measurements are normally made at 3*’hour intervals and any pass 
entering the search area during a period when the buoy was operating 
would do so within 90 minutes of a buoy significant waveheight meas- 
urement. The 90-minute measurement window is Justifiable since wave 
conditions vary slowly. In fact, most of the buoys measured other 
data, such as wind speed, every hour. However, NOAA determined 
that, due to power constraints, it was sufficient to measure SWH 
every three hours (Steele, 1980). It was assumed, however, that no 
storm fronts entered the search area in the interval between the 
time of the buoy measurement and the time of the altimeter estimate. 

After the altimeter passes which entered one or more of 
the search areas within 90 minutes of a buoy measurement were iden- 
tified, it was necessary to ascertain the status of the altimeter 
data. When the altimeter was not locked in the tracking mode (this 
Usually occurred during and shortly after the time when the satel- 
lite passed over land), the data segment could not be Included in 
the set of data to be compared with the buoy data. Had any of the 
SWH estimation algorithms which process IRS data been used, the 
altimeter data would also have had to be in TM Mode 2 with all IRS's 
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reported. This would have reduced the size of the data set by more 
than half. 


The 5WH estimate which occurred at the point of closest 
approach to the buoy was used in tne comparison study. As discussed 
in Chapter 2 , the GEOS-3 significant waveheight estimates were based 
upon a 21-second sliding rectangular filter of the ARS's and the 
ARS's Were averages of the instantaneous waveforms, which were 
reported 100 times per Second. Therefore, the SWH used for compari- 
son with the buoy mea.^urements Was based upon over 2100 instantane- 
ous waveforms. T>/is smoothing removed most of the variability in 
the estimate in the vicinity of the point of closest approach (PCA) 
to the buoy. 

The same mode and status requirements were employed in 
identifying those segments of altimeter data to be included in the 
global atlas of significant w/aveheight. The proximity requirement 
in time and location to a NOAA buoy, however, did not apply to the 
part of the investigation concerned with generating the global atlas 
of significant waveheight. Therefore, the size of the data set used 
for the global atlas was vastly larger than the data set used for 
the buoy comparison. The ^ize of the buoy comparison data set was 
126 passes, while the size of the data set used for the global atlas 
was approximately 8000 passes. The actual processing, which 
included several other algorithms besides the SWH estimation algo- 
rithm, required the use of two computers, operating 2^ hours per 
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day, 7 days per week, for several months. 


3.2 Buoy Data Set 

During the preliminary stages of the investigation, a can- 
puter program was written to select all GEOS-3 altimeter passes 
whose ground tracks passed near one of the NOAA buoys (see Figure 
1.4). ^>e data from those buoys (Hadsell, 1974) which matched the 

GEOS-3 altimeter passes was then identified and requested from NOAA. 
Unfortunately, at this early stage, a search area radius of 1/4 
degree was used to identify the desired buoy data. Some time later, 
it was determined that a one-degree search area could have been used 
to select the buoy data. Although the search area was changed to 
one degree and additional matches were found to have occurred 
between the buoy data already received and the entire GEOS-3 altime- 
ter data set, the number of matching passes of altimeter data and 
buoy data would have been increased by approximately 300 t if the 
requested buoy data had been based upon a one-degree search area 
ladius. However, due to budgetary constraints, it was not possible 
to obtain additional buoy data at a later time. 

The NDBO measurements of significant waveheight were made 
by three different types of buoys. These types were: 

1. EEP - Engineering Experimental Phase 
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2, WSA - Wave Spectrum Analyzer 

3. WDA - Wave Data Analyzer 

The Engineering Experiment Phase type buoy was the first system 
developed by NDBO and is described in Steele, et al (1975). 
Although the system performed acceptably, the hardware was not 
operationally reliable, and the measurements displayed relatively 
high noise levels (Steele and Johnson, 1977). Currently only one 
EEP type buoy is deployed, and that buoy is used only for the pur- 
pose of evaluating the performance of the more advanced WSA and WDA 
type buoys. Only one EEP measurement was used in this investiga- 
tion , 


Most of the buoy data used in this Investigation (about 
80%) were reported by WSA type buoys, which are described in Remond 
(1976). WSA buoys consist of 12 analog filters, each with a 
separate center frequency. The output from a strapped down 
accelerometer is fed into the filters, and a system of 12 equations 
with 12 unknowns is solved to produce the estimate of r ^nificant 
waveheight. However, the filters were not precis' ly calibrated and 
tended to estimate significant waveheight which was 15% 1 wer than 
the true waveheight (Steele, 1980). The only other measurement 
error thought to be of consequence <s the measurement error of the 
accelerometer (Steele, 1980). NDBO estimates that the WSA determi- 
nation of significant waveheight has a measurement error of 50 to 
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100 cm (Withee and G\assingame, 1976). 

About 20} of the buoy data used in this investigation were 
reported by WDA type buoys, which are currently the most advanced 
buoys deployed by NDBO. These buoys are described by Steele, et al 
<1976). In this system, data from a vertical strapped do«vn 
accelerometer is passed trirough an annlog filter, digitized and 
transformed into the equivalent of covariances (Steele, et al, 1975, 
and Steele, et al , 1976), from which significant waveheight is cal- 
culated. NDBO is very Cv‘»ni‘ident in the quality of the WDA measure- 
ments, but Steele (198C) estimates that the significant waveheight 
reported by the WDA type buoys is 10 to 20 cm low. ..s was the case 
with the WSA type buoys, the only other significant error source is 
thought to be the measurement error of the accelerometer (Steele, 
1980). NDBO estimates that the WDA determination of significant 
waveheight has a measurement accuracy of 30 to 50 cm (Steele, 1980). 

Many of the altimeter passes which entered one of the 
search areas could not be matched to buoy data. This occurred 
because the buoys were often deployed for several months, then 
recovered and not deployed again until some time later (see Figure 
3.2). Additionally, more than half of the buoy data was reported 
containing atmospheric measurements, such as wind speed, temperature 
and pressure, but not significant waveheight measurements. Tibbie 
3.1 lists all of the passes which satisfy both the altimeter data 
set requirements and which were requested and received from the NOAA 
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Data Buoy Offict. 


It should be noted that NOAA’s philosophy is (Steele, 
1980) that a high degree of accuracy in measuring sea state is not 
as important for low sea state as it is for high sea state. For 
example, an estimate of SWH = 2 meters in error by a meter or more 
is less important than an estimate of SWH = 5 meters which is in 
error by a meter or more. Additionally, ships and ocean platforms 
are designed to be able to endure certain maximum stresses. Since 
high sea state may cause stresses that approach structural limits, a 
precise determination of SWH in high sea state conditions is impera- 
tive. 


3. 3 Seasat Data Set 

During September and October of 1978, the Gulf of Alaska 
Seasat Experiment (GOASEX) was performed to aid in the accuracy 
determination of the Seasat geophysical parameters. Part of that 
experiment contained altimeter parameter accuracy determination. 
The results of that determination were published in the Seasat Gulf 
of Alask ^ Workshop II Report (1979). In that report, 17 ground 
track crossings of Seasat and GEOS-3 "were compared for the purpose 
of verifying the accuracy of the significant waveheight determina- 
tions of each of the altimeters. The Seasat and GEOS-3 data for 
those 17 passes make up the Seasat data set, which is examined in 
Chapter M. 
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CHAPTER ^ 


COMPARISON OF GEOS-3 SWH ESTIMATES WITH BUOY AND 
SEASAT SWH MEASUREMENTS 

The GEOS-3 significant waveheight estimates and the NOAA 
buoy significant waveheight measurements constitute two independent 
determinations of the same phenomenon. Some information about the 
accuracy of the buoy estimates is available from NOAA, and these 
data, combined with statistical analysis of the difference between 
the GEOS-3 and the buoy determinations, can be used to infer the 
accuracy of the GEOS-3 significant waveheight estimate. 

The analysis that follows assumes that the two estimates 
are independent and that they have independent error sources. 
Although there is no reason to believe that these assumptions are 
violated, it should be noted that the statistical analysis which 
determinr(?s the accuracy of the GEOS-3 significant waveheight phenomenom 
could be corrupted slightly if the GEOS-3 and buoy determinations 
had error contributions from the same indirect source. An example 
of this phenomenom is the degradation of the buoy measurement due to 
wind speed and direction. Since the wind affects the motion of the 
buoys and since the significant waveheight measured by the buoys is 
related to the motion of the platform, both wind speed and direction 
could cause an error in the measurement of significant waveheight by 
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the buoy. Although the GEOS-3 estimate was assumed to be Indepen- 
dent of wind effects, the existence of such dependence, along with a 
similar dependence by the buoy data, would influence the altimeter 
SWH accuracy estimates which follow. 

Table 4.1 lists all of the altimeter passes which could be 
matched with buoy data. If the standard deviation of the differ- 
ences betwet^n the GEOS-3 significant waveheight estimates and the 
NOAA buoy significant waveheight measurement.” is computed, it can be 
seen that the difference between the two estimates for revolution 
number 10227 is noticeably larger in absolute value and in fact is 
greater than three times the standard deviation (3-sigma) of the 
differences of the entire comparison data set. It can therefore be 
legitimately edited from the data set leaving 125 passes of data on 
which the following analysis is based. 

4 . 1 Computation of the GEOS -3 SWH Stands) d Deviation 

As has been stated earlier, one of the primary objectives 
of this investigation is to establish the accuracy of the Wallops 
SWH estimation algorithm. Table 4.1 presents the differences 
between the SWH computed using that algorithm and the SWH measured 
by buoys. Using the statistics of those diferences and the statis- 
tics of the buoy measurement permits the statistics of the GEOS-3 
estimate to be inferred. 
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0]UG!t!AL page IS 
or POOR QUAUTY 


CASE 

REV 

UNIO 

BUOY 

DATE 

1 

681 

626 

E013 

760521 

2 

1111 

86 

EB03 

760827 

3 

1438 

346 

EB03 

760720 

4 

2291 

626 

EB03 

760018 

5 

2380 

62 

EB16 

760026 

6 

2870 

M2 

EB41 

761020 

7 

3032 

603 

EB41 

761110 

8 

3060 

617 

EB41 

761112 

0 

3M1 

780 

EB03 

761203 

10 

3473 

841 

EB16 

761211 

It 

3668 

U2 

EB41 

761217 

12 

3506 

M5 

EB41 

761220 

13 

4213 

423 

EB04 

760201 

14 

4626 

200 

EB16 

760223 

16 

4500 

240 

EB16 

1^228 

18 

4017 

127 

EB16 

760322 

17 

4021 

132 

EB03 

760322 

18 

6116 

MO 

EB16 

780405 

10 

6138 

414 

EB41 

760407 

20 

6443 

227 

EB16 

760428 

21 

6028 

121 

EB41 

760808 

22 

6367 

166 

EB01 

760703 

23 

6616 

120 

EB41 

760713 

24 

6562 

160 

EB41 

760716 

26 

6623 

120 

EB16 

760721 

26 

6629 

135 

EB16 

760721 

27 

6870 

186 

EB41 

760808 

28 

7021 

127 

EB16 

788818 

20 

7867 

176 

EB70 

761017 

») 

7874 

179 

EB16 

761017 

31 

7808 

118 

EB03 

781010 

32 

7024 

146 

EB70 

761021 

33 

7036 

154 

EB78 

761021 

34 

7038 

160 

EB43 

761022 

36 

7038 

160 

EB70 

761022 

38 

7087 

187 

EB10 

761024 

37 

6008 

127 

EB16 

761027 

38 

8066 

181 

EB70 

761031 

30 

8270 

102 

EB71 

761116 

40 

8422 

136 

EB03 

761126 

41 

8606 

126 

EB71 

761208 

42 

8840 

160 

EB70 

761211 

43 

8660 

168 

EB70 

761212 

44 

8848 

166 

EB70 

761226 

46 

8860 

162 

EB70 

761226 

46 

8861 

165 

E616 

761226 

47 

8877 

107 

E^IO 

761227 

48 

8046 

168 

EB41 

770101 

40 

6048 

160 

EB21 

770101 

60 

0076 

106 

EB04 

770110 

61 

0076 

108 

EB03 

770110 

62 

0120 

136 

EB03 

770114 

63 

0132 

138 

EB71 

770114 


PCA 



8WN (Ml 


TIME 

KM 

QEOS 

BUOY 

OIFP 

13720 

64,0 

.M 

1,06 

-.56 

12,1453 

30.0 


.00 

-.87 

161. T 20 

.6 

.00 

,00 

-.00 

222(33 

28,1 

2.21 

2.a 

.01 

204060 

11.1 

3.06 

2.18 

.a 

204030 

a, 2 

1.» 

1,10 

,a 

70332 

13.1 

1,87 

,47 

i,a 

221652 

21.4 

i.a 

1.10 

.a 

131667 

3.0 

2.04 

3.16 

-.22 

110818 

a, 2 

2.11 

4.13 

-2.02 

112200 

14,4 

2.06 

1.10 

,H 

23431 

22.0 

2.a 

1 ,W 

,M 

1826 M 

22.1 

6,60 

4,37 

i.a 

104520 

10.3 

4,45 

3,17 

i.a 

103154 

a.o 

,00 

,70 

-.70 

131613 

a, 2 

1,80 

i,a 

,63 

105351 

21.2 

3.27 

3.02 

a 

146020 

27.3 

4.54 

2.40 

2.06 

1736 

a .6 

2,14 

.60 

1,45 

173336 

22,3 

1,62 

.60 

1,02 

223120 

1.7 

2.56 

.60 

1 ,W 

10066 

12,1 

1,31 

1.83 

-.62 

113720 

6,6 

1,31 

1,73 

-.42 

24041 

12.7 

,44 

.60 

-.06 

31814 

105.0 

1.23 

,61 

.62 

126836 

18,8 

1.77 

,67 

1.10 

63240 

a.o 

1.77 

1,00 

,60 

62816 

11.6 

1,00 

1 ,M 

.10 

11634 

6,3 

3,60 

3.67 

-.17 

132027 

21,0 

1,46 

.80 

.66 

22040 

0,0 

2.01 

1,40 

1,42 

16746 

107,1 

1,66 

1.87 

-.21 

204806 

10,3 

3,61 

1.60 

1,02 

14326 

a ,8 

2.66 

1,60 

,87 

14324 

63,0 

2.66 

1.64 

,02 

26100 

a.o 

2.01 

1,81 

1.10 

1766 

2,0 

2.74 

2,08 

.66 

26130 

a.B 

8 ,a 

8,06 

.12 

35012 

23,1 

2.a 

1,00 

.a 

64441 

17,9 

4.33 

2.03 

1,40 

64213 

33.0 

1,87 

1 .M 

-,01 

76046 

76,4 

2.72 

3.08 

_,38 

24102 

10,4 

6.30 

6,63 

-,M 

02642 

62.1 


4.50 

1.27 

41650 

33.6 

-OBk 

3 .M 

-.14 

71M4 

a ,6 

3.17 

2.04 

1.13 

10 M 37 

a .6 

3.63 

3.76 

-.12 

73217 

a .3 

4,03 

2,86 

1,18 

106817 

27.3 

2.71 

1,60 

1,11 

101820 

a.8 

4.W 

3.60 

.61 

121042 

27,2 

3.66 

3 .a 

.37 

61850 

a.i 

3.86 

3.66 

.21 

110007 

16.7 

2.a 

1.63 

,72 


TABLE 4.1, QEOS/BUOY COMPARISON RESULTS 
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OF FOOR QUALITY 


CASE 

REV 

64 

0188 

65 

0218 

68 

0253 

67 

0273 

58 

0310 

69 

0366 

60 

0044 

61 

0700 

62 

9838 

63 

9008 

64 

10030 

66 

10067 

66 

10080 

67 

10110 

68 

10110 

60 

10128 

70 

10128 

71 

10167 

72 

10163 

73 

10166 

74 

10167 

75 

10208 

76 

10227 

77 

10240 

78 

10266 

70 

10306 

80 

10334 

81 

10360 

82 

10M0 

83 

10384 

84 

10433 

86 

10430 

86 

10633 

87 

10636 

88 

10507 

80 

10632 

00 

10664 

01 

10664 

02 

10703 

03 

10820 

94 

10831 

06 

10060 

96 

10066 

97 

10f179 

98 

11367 

90 

11606 

100 

11642 

101 

11586 

102 

11633 

103 

11741 

104 

11632 

106 

11886 

106 

11883 


UNIQ 

BUOY 

114 

EB70 

141 

EB19 

167 

EB16 

108 

EB41 

134 

EB41 

170 

EB16 

174 

EB43 

138 

EB41 

187 

EB41 

132 

EB41 

172 

EB70 

187 

EB16 

116 

EB16 

140 

EB03 

140 

EB17 

153 

EB03 

163 

EB16 

174 

ES17 

103 

EB41 

106 

EB10 

\Of, 

EB70 

134 

EB16 

148 

EB70 

167 

EB16 

166 

EB16 

134 

EB16 

161 

EB71 

106 

EB43 

106 

EB70 

117 

EB16 

146 

EB15 

140 

EB16 

160 

EB71 

163 

EB16 

120 

EB19 

166 

EB16 

101 

E803 

101 

EB16 

134 

EB01 

147 

EB10 

163 

EB16 

108 

EB16 

112 

EB16 

123 

EB01 

107 

EB16 

136 

EB41 

164 

EB41 

106 

EB01 

148 

EB01 

163 

EB41 

140 

EB01 

110 

EB04 

130 

ESDI 


PCA 

DATE TIME 

7701 IB 68869 
770120 1W204 
770123 4828 

770124 101539 
7:W127 1 2758 
770130 54858 
770210 154510 
770302 143330 
770305 54M8 
770310 1608^ 
770310 136045 
770320 201440 
770322 113940 
770324 142750 
770324 142032 
7703» 204623 
770325 204131 
770327 216637 
770328 82843 
7703» 13»48 
770328 160747 
770331 124742 
770401 204421 
770402 183733 
770403 214033 
770407 02404 
770400 103742 
770411 104234 
770411 164238 
770413 226734 
770416 103206 
770416 2012K 
770423 121233 
770423 163035 
770428 2034 

770430 120666 
770602 10406 
770602 6014 

770506 123241 
770613 166631 
770614 134146 
770623 144046 
770624 3006 

770626 1666 

770620 175021 
770701 43433 
770703 104660 
770704 193306 
770710 64120 
770717 212136 
770724 71614 

770726 224863 
770727 221640 


KM 

OEO) 

15.6 

5.21 

30,1 

4.60 

33,0 

4.38 

28.3 

.87 

17,1 

3,36 

20.1 

4.40 

7,7 

6.18 

0.8 

1.82 

.0 

3.23 

45.4 

1.80 

6.8 

1.66 

44.7 

1,14 

84.6 

2.10 

36,7 

2.87 

12.8 

2.65 

8.3 

3.00 

02.6 

2.18 

46.7 

Z45 

30,1 

.00 

6.7 

4.69 

36.2 

2.23 

2.4 

1.46 

66.0 

,00 

37.0 

2.04 

9.3 

2.31 

6.0 

2.66 

60.7 

2.26 

12.6 

2.07 

10,8 

2.07 

77,6 

6.7R 

90,8 

2.®1 

24.6 

2.64 

14.6 

.00 

36.P 

3.02 

27,0 

1,13 

38.2 

1,06 

2.6 

3.33 

109.7 

1.40 

33.7 

1.17 

46.4 

2.69 

24.6 

1,13 

78.7 

3.66 

46.7 

3.66 

50,1 

2.26 

46.0 

1.80 

46,4 

1.74 

36.3 

,00 

17.3 

1.18 

27.6 

2.4 6 

20.6 

.40 

3bA 

1.14 

12,0 

1.10 

21.4 

2.16 


VNH |M| 


BUOY 

DIPP 

5.07 

-.78 

3.87 

.73 

3.23 

1,16 

1.76 

.-.88 

2.56 

.79 

3.n 

.42 

8.77 

.41 

1.80 

.02 

3.12 

,11 

.64 

1.05 

1.37 

.20 

1.53 

-.30 

1.88 

.31 

1.03 

.04 

2.60 

-.06 

3.02 

.07 

1.72 

.48 

2.73 

-.a 

.70 

-.70 

3.06 

,64 

2.10 

.04 

1,00 

-.44 

2.36 

-2.» 

1.64 

.40 

1.80 

.61 

1.17 

1.30 

1.10 

1.18 

2.83 

.14 

3.00 

-.03 

3.87 

1,88 

1.00 

1.81 

1.13 

1.61 

.81 

-.81 

2.83 

1,00 

1.») 

-.17 

1.30 

-.24 

3.10 

.14 

1.13 

.27 

1,34 

-.17 

3.10 

-.60 

1.11 

.02 

2.20 

i.a 

2.62 

.93 

1.40 

.76 

1.40 

.40 

.06 

.70 

,84 

-.84 

1,06 

.13 

1.20 

I.a 

1.00 

-.61 

1.72 

-.68 

,60 

.50 

1.68 

.47 
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' . . . I I ilab, 

.Of- POOR QUALITY 


CASE 

REV 

UNIO 

BUOY 

DATE 

RCA 

TIIME 

KMI 

OEM 

SWHOMI 

BUOY 

Dirr 

10? 

12028 

171 

EB20 

77M07 

41182 

44.4 

1.88 

M 

.73 

10B 

12031 

174 

EB41 

770807 

W|42 

17.8 

1.67 

1.«> 

.87 

100 

12M8 

121 

EB41 

770010 

417 

13.8 

1,27 

1.00 

,27 

110 

12188 

121 

EB01 

770818 

98838 

7.8 

237 

1.03 

1.34 

111 

12210 

177 

EB18 

770820 

8818 

38.8 

288 

1.S) 

1.88 

112 

12217 

184 

EB71 

770820 

121»1 

48.8 

1.80 

.80 

1.30 

113 

1U30 

188 

EB41 

770821 

102828 

42.3 

.81 

1.00 

-.38 

114 

12 »e 

180 

EB41 

770802 

24888 

48.0 

vOO 

.80 

-.80 

lie 

12408 

182 

EB01 

770803 

23314 

8.7 

1.87 

.88 

.88 

lie 

124M 

141 

EB01 

770808 

124137 

37.2 

204 

1.00 

1.04 

117 

12871 

138 

EB18 

770821 

143028 

21.8 

218 

,82 

1.28 

118 

12084 

128 

EB18 

771008 

82482 

38.2 

207 

1.37 

.70 

lie 

12888 

134 

EB17 

771008 

130812 

3.0 

1.08 

1.43 

-.38 

120 

12870 

137 

E»18 

771008 

180812 

WA 

.81 

.83 

-.02 

121 

12tt4 

182 

EB41 

771008 

188208 

30.8 

1.82 

.88 

.78 

122 

12821 

188 

EB41 

771008 

7042Y 

20.8 

1.88 

1.44 

.48 

123 

12888 

187 

EB18 

771014 

171301 

13.7 

230 

270 

.80 

124 

13083 

128 

EB18 

771018 

78833 

27.3 

1.83 

1.03 

.80 

12B 

13120 

182 

EB41 

771023 

83802 

38.1 

.00 

.83 

-.83 

120 

13187 

181 

EB16 

771028 

184747 

80.8 

.00 

1.17 

-1.17 
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Using an edit criterion of three times the Standard devia- 
tion of the differences between the altimeter and buoy determina- 
tions of significant waveheight, the standard deviation of the 
differences not edited (125 passes) is 

Sp = 0.71 m (4.1) 

The letter S is used to symbolize standard deviations of the sam- 
ple population of 125 elements and must be differentiated from the 
symbol a , which will be used to symbolize the standard deviation 
of the entire population. 


With the assumption that the GEOS-3 estimates and the buoy 
measurements of significant waveheight are independent, then 


a 


2 

D 



(4.2) 


where cTg is the standard deviation of the buoy significant 
Waveheight measurements, a„ is the standard deviation of the altiro- 

U 

eter significant waveheight estimates and is the standard devia- 
tion of the difference between the GEOS-3 estimate and the buoy 
measurement of significant waveheight. The value of is an 

estimate of and can be substituted into Equation (4.2) to yield 


S 


2 

D 



(4.3) 


Equation (4.3) could be used to solve for if Og was known. 

The NOAA Data Buoy Office estimates that 0.50 £ Og £ 1.00 m 
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(Withee and Blassingame, 1976). 


Examination of Equation (4.3) shows that as 


a-g 1. 

decreases. Therefore, the 

is 

used. Substituting this value and Equation (4.1) into Equation 
(4.3) yields 


increased from 0.50 m to 1.00 m, 0 . 


most conservative estimate for Oq is obtained when Og - 0.50 


Sq = 0.50 m 


(4.4) 


where the symbol S is again used to indicate that the standard 
deviation refers to the sample population. 

The standard deviation of the entire population can be 
estimated to the 95f confidence limits by employing the chi-square 
distribution. These limits are given by 




(4.5) 


where n is the number of degrees of freedom (in this case 

2 

n s 123). When n > 30 , a normal approximation can be used and 
can be determined from 


= I (X(j ♦ [2n-1]^^^)^ (4.6) 

where is the a-point of the cumulative normal distribution. 

From the normal probability function tables 
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F(ic) s 0.025 


(4,7) 


*,025 ^ 


F<y.) = 0.975 


*.975 ■ 


(4.8) 


Substituting Equations (4.7) and (4,6) into Equation (4.6) and using 
n s 123 yields 


^.025 




.975 


= 93.7 


= 155.1 


(4.9) 

(4.10) 


Substituting Equations (4,9) and (4.10) into Equation (4.5) yields 


0.33 >^Q> 0.20 m 


(4.11) 


and since cIq must be non-negative 

0,57 > > 0.45 m (4.12) 

ij — 

Therefore, based upon the assumption of Independence between the 
altimeter estimates and the buoy measurements of significant 
waveheights, the 3-slgma edit crlteron, a value of cJg = 0.50 m , 
and the sample population of 125 differences, there is a 95t proba- 
bility that the value of for the entire population of' GEOS-3 

passes lies between 0.45 m and 0.57 m. 


The above values represent a conservative estimate under 
the assumptions given above. To Illustrate, it can be shown that 
employing a 2 1/2-sigma edit criterion results in an estimate of 
0,40 m for <y„ with a 95% probability that the standard deviation 

(j 
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of the entire population of GEOS-3 passes lies between O .36 m and 


0,46 m. Similarly, the estimate of 0 ^ would be lowered if the 

estimate of 0 was chosen to be larger than the minimum estimate 
21 

published by NOAA. Based upon these results, the GEOS-3 design 
specification of+ 20 % accuracy for 2.0 < SWH < 10 meters appears to 
be satisfied. 

Table 4.2 lists the statistics of the comparisons of 
GbOS-3 estimates and buoy measurements for each of the buoys. Note 
that the results obtained for buoy EB15 are significantly worse than 
the results obtained for the other buoys. Due to this result, NOAA 
is currently investigating the accuracy of the significant 
waveheight measurf'inents of buoy EB15. 

4.2 Computation of the SWH Mean Difference 

Using the same 3-*sigma edit criterion for the difference 
between the GEOS-3 significant waveheight measurements that was 
described in Section 4.1, the mean difference of the 125 samples is 

D = 0.41 m (4.13) 

In Equation (4.13), the letter D is used to symbolize that the 
mean difference given is for the sample population. The mean 
difference for the entire population will be given the symbol . 
Figure 4.1 Illustrates that the differences are distributed normally 
about the mean of 0.41 m. It should be noted that the NOAA Data 
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BUOY BUOY TYPE MEAN SIGMA # SAMPLES 
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FIGURE 4 . 1 . OISTRimJTlOW OF SWH, 


Buoy Office (Steele, 1980) estimates that they probably over- 
corrected for noise, especially for low sea state and that the over- 
correction would tend to make the buoys measure 15% low. This is 
consistent with the mean difference given in Equation (<4.13). 


The mean difference of the entire population can be com- 
puted to the 95% confidence limits by employing the t-distr ibution. 
These limits are 


D - 


^. 975^0 


Cn] 


1/2 - 




D + 


^. 975^0 


Cn] 


1/2 


(4.14) 


From the standard t-distr ibution tables, 

t 975 = ■)-980 (4.15) 

for n = 123. Substituting Equations (4.1), (4.13) and (4.15) into 
Equation (4.14) yields 

0.32 < Vp < 0.50 m (4. 16) 

Therefore, based upon the assumption of Independence between the 
altimeter estimates and the buoy measurements of significant 
waveheight, the 3-sigraa edit criterion and the sample population of 
125 differences, there is a 95% prob hility that the value of Pp 
for the entire population of GEOS-3 passes lies between 0.32 m and 
0.50 m. 


As was the case vHh o„ , the value of the mean differ 

U 
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ence for the population can be made to decrease by changing the edit 

criterion. It is not, however, a function of the estimated buoy 

standard deviation o . Although the magnitude of the mean differ- 

u 

ence cannot be precisely determined from this data set, that any 
mean difference is present is, in itself, an important result. 
Furthermore, as was pointed out earlier in this section, the mean 
difference could be entirely accounted for by the NDBO overoorrec- 
tion for noise. 

Two points should be cited here in connection with the 
mean difference. First, the altimeter significant waveheight esti- 
mate was calibrated by setting the smallest expected value of the 
•lignificant waveheight algorithm to zero. This calibration is sub- 
ject to error and could easily account for a small bias between the 
two estimates. 

Second, the data which account for most of the mean 
difference are the data where significant waveheight is small. This 
is due to the previou.«£ly mentioned algebraic sensitivity of the 
altimeter estimation algorithm for near-calm seas and to the fact 
that the magnitude of low sea states is often smaller than the pre- 
cision of the buoys. This is readily seen by examining Table 4.1. 

4 . 3 Linear Regression 

Since the GEOS-3 significant waveheight estimate jjijo “ft*? 
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NOAA buoy significant waveheight measurement constitute two indepen- 
dent determinations of the same observable quantity, a linear 
regression analysis may be performed upon the two sets of data to 
determine their agreement. If both the altimeter estimate and the 
buoy measurement were perfect, a positive correlation with unity 
magnitude would be expected. The degree to which the actual corre- 
lation coefficient differs from unity is a measure of relative 
errors of the two estimates. 


For a set of paired data points (x^^, y^, i n t , 2, 3, • . . ,N) , 
the line which best fits the data (in the sense of least squares) is 
given 

y=mx+b (M.17) 

where m is the slope of the line and b is the y-intercept of the 
line with m and b given by 


m 


N 

1 = 1 


NO^ 


- X y 


('<.18) 


b = y - mx 


(4.19) 


where x and y are the mean values of x^ and y^ , respec 
tively. The correlation coefficient X is given by 


y = 


ma 
X 

a 

y 


(4.20) 
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a and 
X 


o 

y 


are given by 
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-2 
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i = 1 



(4.21) 


(4.22) 


Additionally, the standard error of estimate, , is computed from 


Sg. = Oy [1 - 


(4.23) 


Figure 4.2 represents a linear regression analysis of the 
data presented in Table 4.1. The GEOS-3 significant waveheight 
estimate is the dependent variable, and the buoy significant 
waveheight measurement is the independent variable. Using the 
analysis presented above 


m 

= 0.96 

(4.24) 

b 

= 0.50 

(4.25) 

y 

= 0.86 

(4.26) 


r 0.75 

(4.27) 


Similarly, Figure 4,3 represents a linear regression 
analysis of the data but with the altimeter significant waveheight 
vistimate as the independent variable and the buoy significant 
waveheight measurement as the dependent variable. Under these con- 
ditions. 
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m 

= 0.77 

(4, 

.28) 

b 

= 0.13 

(4- 

.29) 

y 

s 0.86 

(4. 

,30) 

Sg 

= 0.67 

(4. 

,31) 


The correlation coefficient given in 

Equations (4.26) 

and 


(4.30) indicate that the two data sets are significaf«tly correlated, 
as would be expected. Although it is difficult to make any more 
than a qualitative assessment from the two linear regressions, it is 
important to remember that regression analysis assumes that the 
independent variable is without error. With this in mind and using 
the fact that the slope in Figure 4.2 is very close to unity and the 
slope in Figure 4.3 is slightly degraded, it might be inferred that 
the errors in the GEOS-3 significant waveheight estimates are 
slightly more prominent than the errors in the buoy significant 
waveheight measurements. Nevertheless, the important result of this 
analysis is the high correlation between the two data sets. 

Figure 4.4 illustrates the difference between the altime- 
ter estimate of significant waveheight and the buoy measurement of 
significant waveheight as a function of the distance from the point 
of closest approach to the buoy. Because the original data request 
made to the NOAA Data Buoy Office was based upon a 1/4 degree search 
area radius, there are more points with small PCA distance than with 
large PCA distance. Nevertheless, Figure 4.4 shows that the differ- 
ence between the two estimates of significant waveheight is not a 


85 



86 


ruhction or PCA dlatance up to at least one equatorial degree (111 
kilometers) . 

GOASEX Comparison 

The Guir of Alaska Seasat Experiment (GOASEX) was designed 
to aid in the determination of the accuracy of the Seasat goophysi> 
cal parameters, including significant waveheight ( Seasat Gulf of 
Alaska Workshop II Report , 1979). In that report, the significant 
wave-heights from 17 ground track crossings of GEOS-3 and Seasat were 
compared and their statistics presented. For those 17 passes, the 
report found that the me.jin difference between the GEOS-3 real-time 
and the Seasat on-board estimates of significant waveheight was 59 
cm. The GEOS-3 data used for that comparison was obtained from the 
GEOS-3 Near-Real-Time Data System (McMillan, 1978). 

It should be noted that the GEOS-3 and Seasat data used in 
the comparison were not published in the report, but that only the 
statistics of that comparison were. For the purposes of verifying 
the results of that comparison, the data set was obtained from S. 
Fedor (1981) of the NOAA Environmental Research Laboratory, who 
served as chairman of both the GOASEX Altimeter Wind-Wave Panel and 
the GOASEX Sensor Intercompar ison Panel. Examination of the data 
employed in the comparison led to the conclusion that the GEOS-3 
data used had been Improperly selected. The GEOS-3 project scien- 
tist, H. Ray Stanley (1981), concurred with that conclusion. The 
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proper GEOS-3 comparison data set is given in Table 

Analysis of the data in Table ^.3 yields a mean difference 
between the GEOS-3 real-time and the Seasat on-board determinations 
of significant waveheight of 2^ cm, compared to a computed value of 
59 cm published in the GOASEX report. (Note that the mean differ-^ 
ence between the GEOS-3 estimate of significant waveheight and the 
Fedor’ Seasat estimate of significant waveheight is only 15 cm, com- 
pared to 51 cm published in the GOASEX report.) That report con- 
cluded that, 


The important point to note here is that all of the algorithms 
and the Seasat data seem to produce a lower value than the 
GEOS-3 estimates. This is consistent with other observations 
that the GEOS-3 waveheights tend to be slightly biased for the 
low waveheights encountered during this comparison, i.e., 
1 m < H 1/3 < ^ m. 


The above statement Implies that the GEOS-3 determination of signi- 
ficant waveheight was always smaller than the Seasat determination. 
As can be seen from Table ^.3, this is not the case. 


Although the data set presented in Table 4.3 does indicate 
a mean bias between the GEOS-3 and Seasat significant waveheight 
estimates, the magnitude of that bias is considerably smaller than 
the bias reported in the GOASEX report. In fact, the bias is suffi- 
ciently small to be considered a confirmation of the accuracy of 
both determinations. 
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Table 4.3 6E0S-3/SEASAT Comparison Statistics for Sub-Track Crossings 



CHAPTER 5 

GLOBAL ATLAS OF SWH DATA 
5, 1 Description of the SWH Global Atlas 

During the CEOS-3 mission, altimeter data spanning It^OO 
hours were received and processed at Wallops Flight Center. Approx- 
imately 93f of the received data was intensive mode and, therefore, 
suitable for the calculation of significant waveheight. This large 
data set, representing approximately 8,000 GEOS-3 dat:ii segments, was 
used in its entirety in the generation of the significant waveheight 
global atlitS which follows. All 8,000 intensive mode GEOS-3 passes 
were processed by a computer program which used the significant 
waveheight estimation algorithm and smoothing technique presented in 
Chapter 2. That program ran continuously for several months in 
order to accomplish this task. (It should be noted that less than 
40f of the GEOS-3 data set was high data rate intensive mode data 
and, therefore, suitable for processing by the algorithms that use 
IRS data.) 


Due to the vast size of the data set, it was necessary to 
partition the oceans into subdivisions where the statistics of the 
estimated significant waveheight could be accumulated, rather than 
the individual data points themselves. The size chosen for the sub- 
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divisions was a square. TTits was the smallest feasible parti- 

tion size consistent with th? storage capabilities of the computer 
software and the contouring limitations of the plot hardware. 

The significant waveheight statistics for each of the 
1® X 1® partitions were retained for a range of values: 

1. SWH < 1.5 m 

2. SWH < 2.5 m 

3. SWH > 3.5 m 

SWH > 6.0 m 

These partition levels were suggested by the U. S. Navy Marine 
Climatic Atlas of the World (C. S. Navy, 197M). By using the same 
partition levels that were used in that publication, the results of 
the GEOS-3 estimates of significant waveheight could be compared 
with yet another independent source of significant waveheight data. 

Upon examination of the statistics for the 1°x1° subdivi- 
sions, it was found that the number of data points per square varied 
widely from a high of several thousand in the GEOS-3 calibration 
area (bounded by Wali.ops Island, Bermuda, Merritt Island and Grand 
Turk) to a low of seve^^^l dozen in the extreme latitudes of the 
southern Indian Ocean. The da yrseness in certain remote 

regions combined with the fi' t that thf/ statistics of a subdi- 
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vision which contained only a few GEOS-3 data passes could be cor- 
rupted by a single pass of data taken during inclement weather con- 
ditions, suggested that the 1^x1^ subdivisions should be smoothed. 
This was accomplished by combining the results of the 1^x1*^ subdivi- 
sions into larger subdivisions. 

Several smoothing sizes were examined, including 3^x3^, 
5®x5°, and 9°x9°. In each of these cases, all the l®x1® sub- 

divisions in the smoothing area were averaged without weighting. 
After inspection of the smoothed results, the 7°x7° smoothing area 
was chosen since it removed the spurious points without eliminating 
the fine detail near coastal zones. Accordingly, the contour levels 
which fo’low for a given latitude and longitude, <j) and X , respec- 
tively, contain data within the 7°x7® square given by 

(j) - 3 1/2 < ^ < ((I + 3 1/2 (5.1) 

X - 3 1/2 < X < X + 3 1/2 (5.2) 

In addition to the four signific/ant waveheight partition 
levels defined above, it was necessary to represent the temporal 
trends of the significant waveheight data. Therefore, the contour 
maps which follow are segmented in time as follows: 

1. December through February 

2. March through May 
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3. June through August 


September through November 

5. Enti'^e GEOS-3 mission (April, 1975, to December, 1978) 

The first four divisions correspond to the four seasons and show 
global trends in the significant waveheight data on an annual basis. 

All contour maps contain Isopleths of 0-10$, 11-30$, 31-51^$, 51-7'/$, 
71-90$ and 91-100$. 

The GEOS-3 significant waveheight global atl& ; is 
presented in Vol. II, Figures A. 1 to A. 20. It must be pointed out that 
the atlas presented here Is "global" ip the sense that data have been 
compiled from nearly all parts of the world. However, there are 
certain areas of the atlas which contain no contour designation. 
These areas are due to gaps In the data or to holes in the telemetry 
coverage for that particular season. Additionally, since the incli- 
nation of GEOS-3 was 115°, no data were obtained outside the region 
-65° < (^ < 465°. 

5*2 Comparison with the Navy Climatic Atlas 

The comparison between the GEOS-3 derived significant 
waveheight contour maps and the IJ. S. Navy Climatic Atlas (U. S. 
Navy, 197A) contour maps is presented here. In particular, the 
differences between the two sources of significant waveheight for 



the northern Atlantic Ocean will be analyzed because both sources 
are derived from a larger data set in that region and, therefore, 
are presumably more accurate in that region. 

Each of the northern Atlantic Ocean contour maps was com- 
pared with the significant waveheight contour maps published in the 
U. S. Navy atlas. That atlas contains contour maps for each month 
rather than each season so that a direct comparison is difficult. 
The analysis was accomplished by examining the U. S. Navy contour 
'.iap for the middle month in the seasonal GEOS-3 contour map and not- 
ing any significant differences between the first and last month and 
the middle month. For example, the GEOS-3 contour maps for March 
through May (Vol. II, Fig. A. 5 through A. 8) were compared with the _U. S. 
Navy Climatic Atlas contour maps for April, and any significant 
differences between the March and April and the April and May con- 
tour maps in the U. S. Navy publication were noted. 

To illustrate the method by which the two atlases were 

compared, examine Figure 5.1 , which is a reproduction of one of 

the contour maps from the U. S. Navy atlas. Unfortunately, the U. 

S. Navy contour maps were drawn in different colors, which were not 
able to be reproduced in this investigation. However, from examina- 
tion of Figure 5.1 in the region lying within latitudes 10° to f?0° 

North and longitudes 60° to 70° West, it is apparent that the U. S. 

Navy atlas estimates that significant waveheight is less than 1.5 
meters 50% to 60H of the time in the month of January. The December 
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and February contour maps In the U, S. Navy atlis are simUar for 
this area. Volunui 11, A.l, however, Indicates that the altimeter esti- 

mates of significant waveheight in this region were less than 1.5 
meters only 10% to 30% of the time. 

5.2.1 December through February Comparison 

The GEOr)-3 significant waveheight contour maps for the 
months of December through February (Vol. II, Fi^j. A.l through A. 4) were 
compared with the corresponding U. S. Navy Climatic Atlas contours 
in the northern Atlantic Ocean, and the results of that comparison 
are as follows: 

1. For significant waveheight less than 1.5 meters, the GEOS-3 
contours show a lower percentage of data in the latitudes less 
than 30° than do the U. S. Navy contours. The GEOS-3 estimates 
in the South Atlantic Bight indicate 10% to 30% of the data is 
less than 1.5 meters, while the Navy data indicates 30% to 50%. 

In the southern portion midway between the South American and 
African continents, the Navy data Indicates 30% of the signifi- 
cant waveheight measurements are less than 1.5 meters, while 
the GEOS-3 data, though sparse, estimates that percentage at 
less than 10%. 

2. For significant waveheight less than 2.5 meters, the Navy atlas 
again estimates a higher percentage of datv'' near North America 
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than do the GKOS-3 estirietes. In the South Atlantic Dighh, the 
GEOS-3 data indioatea 50f to "JOi of the significant waveheight 
estimates are less than 2.S meters. The Navy atlas estimates 
80$ to 90$ of the waves are less than 3.5 meters in that 
region. Tlje Gulf of Mexico and the Caribbean Sea contour maps 
agree satisfactorily. 

3. For the contour maps of significant waveheight greater than 3.5 
meters, the GKOS-3 estimates indicate a higher percentage of 
values above a latitude of 40°. Tlie GEOS-3 estimates indicate 
that a large region of the northern Atlantic Ocean has signifi- 
cant Wvavetieight greater than 3.5 meters for ?0$ to 90$ of the 
time during these months. Tlie Navy atlas estimates that per- 
centage to be 50$ to 60$. Similarly* the GEOS-3 contours indi- 
cate a region in the South Atlantic Bight and another region in 
the Gulf of Mexico where over 10$ of the significant waveheight 
data is greater than 3.5 meters. Ihe Navy atlas shows no such 
area . 

4. The GEOS-3 contours of significant waveheight greater than 6.0 
meters indicate a higher percentage of data in the extreme 
northern latitudes for these months than does the Navy atl/<s. 
The GEOS-3 contours show a large area where SWH exceeds 6.0 
meters 30$ to 50$ of the time. The Navy atlas indicates SWH 
exceeding 6.0 meters only 15$‘ of the time in that region. 
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Summarizing the November through February data, there are 
fewer calm sea state estimates for GE0S>3 than for the Navy data. 
Conversely, there are more high sea state estimates for GEOS-3 than 
appear in the U. S. Navy Cl imatio Atlas . 

5.2.2 March through May Comparison 

The results of the comparison in the northern Atlantic 
Ocean for the months of March through May (Vol. 11, Fig. A. 5 through A. 8) 
are : 

1. For significant waveheight less than 1.5 meters, the GEOS-3 
data showed a much higher percentage of data in the mid north* 
ern latitudes, while finding a significantly smaller percentage 
of data in the southernmost latitudes of the northern Atlantic 
Ocean. The western coast of the African continent showed con- 
siderable discrepancy, with the GEOS-3 data indicating most of 
the data to be within the 10$ to 30$ range, while the Navy 
atlas found 60$, 70$ and even 80$ of the data to be less than 
1,5 meters in some portions of that area. Both sources agreed 
rather well in the Caribbean, but the GEOS-3 data showed the 
percentage of significant waveheight data in the central Gulf 
of Mexico region to be 30$ to 50$ while the Navy atlas indi- 
cated 60$ to 70$. 

2. For significant waveheight less than 2.5 meters, the GEOS-3 
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data in the middle to northernmost latitudes indicated that SWH 
was less than 2.5 meters lOt to 30^ of the time. The Navy 
atli»is estimated that SWH was less than 2.5 meters to 60^ of 
the time in that region. 

3. The GEOS-3 data indicated a large area in the extreme northern 
latitudes where over 50S and even over 70} of the data was 
greater than 3.5 meters. The Navy atlas showed only 30} of the 
data to be greater than 3>5 meters, except for March when 40} 
of the data exceeded 3.5 meters. All other portions of the 
northern Atlantic Ocean agreed rather well. 

4. For significant waveheight greater than 6.0 meters, the Navy 
atlas showed about 5} in the northern extremities whereas the 
GEOS-3 contour mapt indicated a large area of 10} to 30}. No 
significant amounts of data were found to be greater than 6.0 
meters in the southern part of the region for either source. 

Summarizing the March through May data, the GEOS-3 esti- 
mates of significant waveheight indicate a higher percentage of high 
waves and a lower percentage of low waves in the mid to northern 
latitudes of the region. The southern latitudes of the region agree 
rather well. In the Caribbean and the Gulf of Mexico, the GEOS-3 
estimates are considerably less calm than the results from the U. S. 
Navy Climatic Atlas . 
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5.2.3 June through AuRUst Comparison 

I II I — I ■ nil I •' II— — .^W»WP— 


The results of the comparison in the northern Atlantic 
Ocean for the months June through August (Vol. II, Fig. A. 9 to A. 12) are 

For significant waveheight less than 1.5 meters, the Navy atlas 
showed contour levels of 20 % to 50? north of latitude M5°, 
while the GEQS-3 results indicated that percentage was in the 
range of 10? to 30?. Of even more significance was the South 

Atlantic Bight rea,ion. For this area, the GEOS-3 results indi- 
cated 30? to 50? of the data was less than 1.5 meters, while 

the Navy atlas indicated 70? to 80? was less than 1.5 meters. 

The Navy atlas also showed a higher percentage of calm ,eas in 
the Gulf of Mexico. Here, the Navy estimates that 70? to 80? 
of the significant waveheight estimates are less than 1.5 
meters, while the GEOS-3 result.5 indicate that percentage 
should be 50? to 70?. 

2. For significant waveheight less than 2.5 meters, the two 
sources agree rather well in the mid-latitudes. In the extreme 
northern latitudes the GEOS-3 estimates indicate an area where 
30? to 50? of the data is less than 2.5 meters, while the Navy 
contours show that portion of the area contains 70? to 80? of 
the significant waveheight measurements less than 2.5 meters. 
Similarly, in the &5uth Atlantic Bight region, the Navy con- 
tours indicate over 90? of the data is less than 2.5 meters. 
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while the GEOS-3 data estimatea 701 to 901 of the data to be in 
that range. 

3. The GEOS-3 cohtoura and the U. S. Navy oontoura agree well for 
the percentage of data greater than 3.5 meter a and greater than 
6.0 metera during theae montha. 

Summarillng the June through Auguat data, the percentage 
of GEOS-3 estimatea of significant waveheight are lower than the 
percentage of Navy estimates of significant waveheight for low to 
calm seas during these months. The high seas contour levels are 
similar . 

5.2.^ September through November Comparison 

Tile results of the comparison in the northern Atlantic 
Ocean for the months of September through November (Vol. II, Fit’. A. 13 
A. 16) are: 

1. For significant waveheight less than 1.5 meters, the two 
sources compare rather favorably with the exception of the mid 
latitudes and the extreme northern latitudes. In the mid lati- 
tudes, the Navy estimates that 30} to 60t of the data is less 
than 1.5 meters, while the GEOS-3 data Indicates that most of 
the offshore contours in this area are lOH to 30H, In the 
extreme northern latitudes, the GEOS-3 contours indicate less 



than 10% of the data is lesi, than 1.5 meters, while the Navy 
atlas indicates 10% to 30% of the data is less than 1.5 meters. 

2. The GEOS-3 and Navy contour maps for significant waveheigbt 
less than 2.5 meters agree well in most areas. There is, how- 
ever a significant difference in the extreme northern lati- 
tudes. Thie U. S. Navy Climatic A tlas estimates as high as 80% 
of the data in this region less than 2.5 meters, while the 
GEOS-3 estimates indicate only 10% to 30% of the data has 
values of significant wavehelght less than 2.5 meters. 

3. The GEOS-3 contours of significant waveheight indicate that 50% 
to 70% of the data in the extreme northern latitudes was 
estimated to be greater than 3.5 meters. The U. S. Navy cc.^- 
tours in the same region indicate only 30% to 40% of the data 
to be greater than 3.5 meters. The remainder of the Northern 
Atlantic Ocean contours of GEOS-3 data agreed well with the 
Navy contours. Additionally, the contours of significant 
waveheight greater than 6.0 meters agreed well throughout the 
entire area. 

Summarizing the September through November data, the 
GEOS-3 estimates indicate a much lower percentage of high seas in 
the extreme northern latitudes than did the Navy atlas. The GEOS-3 
data also indicated fewer calm seas off the U. S. coast than did the 
y* ^3vy Climatic Atlas . 
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5.3 Discus sion of the Contoujp Compari sons 

From the information presented in Sections 5.2.1 through 
5,2.4, one overriding point is clear: there are significant differ- 

trces between the U. S. Navy Cl im '^tic Atlas and the GEOS-3 signifi- 
cant waveheight atlas. This is an extremely important result which 
could have significant impact on such areas as ship routing and 
design of offshore construction. These GEOS-3 results have poten- 
tial for providing an overall improvement in significant waveheight 
atlases, especially when combined with other sources to provide data 
for the poorly determined GEOS-3 regions. 

The major differences between the two contour atlases can 
be summarized as follows: 

1. The GEOS-3 contours indicate a significantly smaller percentage 
of calm seas in the lower latitudes. 

2. The GEOS-3 contours in the extreme northern latitudes indicate 
a significantly higher percentage of high sea states. 

3. The Navy contours indicate a higher percentage of significant 
waveheight measurements in the range of 2.5 to 3.5 meters. 

It should be pointed out that since the GEOS-3 estimate 
has been proven to be accurate by comparison with buoy and Seasat 
d?nba, the Navy contour maps should at least be scrutinized care- 
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fully. The differerices between the two atlases could arise because 
of the error sources in each (sparse measurements in various 
regions, instrument errors, errors caused by human measurements, 
etc.) or because the GEOS-3 contour estimates are made up of only 
four years of significant waveheight data, whereas the Navy data has 
been compiled over decades, 

C, L. Parsons, a NASA meteorologist at Wallops Flight 
Center who has published several papers on the GEOS-3 significant 
waveheight estimate, believes (1976) that these differences are not 
unexpected. The U. S. Navy data used in the extreme latitudes are 
quite sparse, especially during the winter months. Additionally, 
the Navy atlas, in its preface, states: 

Waveheights reported by most transient ships tend to be 
underestimated in comparison to those recorded by Ocean Weather 
Stations and dictated by the synoptic situation. 

This is consistent with the conclusions stated earlier in this sec- 
tion that the U. S. Navy contour maps indicate a higher percentage 
of low sea states and a lower percentage of high sea states than do 
the GEOS-3 contour maps. 

In Parsons (July, 1979), it is demonstrated that the 
GEOS-3 significant waveheight measurements agree better with wave 
recorder data than do measurements made by a human observer. Par- 
sons (September, 1979) also points out that ships avoid high sea 
state areas where possible and shows storm tracks where GEOS-3 meas- 


urements agreed well with ship measurements, but indicated that very 
Tew ship measurements were available within the most violent areas 
of the storms. 

These factors, combined with the fact that all of the Navy 
data in the extreme latitudes come from human estimates made onboard 
transient ships, suggests that the GEOS-3 estimate would be higher 
in these regions. Parsons also believes (1978) that the GEOS-3 
estimates contain fewer calm seas near the coastal regions of North 
America because of the numerical problems encountered by the altime- 
ter estimation algorithm in regions of calm seas (see Section 2.5). 
H. Ray Stanley, the GEOS-3 project scientist, concurs in both of 
these conclusions. 
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CHAPTER 6 

¥■ 

CONCLUSIONS AND RECOMMENDATIONS 


An algorithm has been derived for computing ocean signifi- 
cant wavehelght from the altimeter measurements of the GEOS-3 satel- 
lite. The al» ' meter-derived estimates appear to be accurate to 
about 50 cm, ^..sod upon a comparison of 125 statistically represen- 
tative passes of GEOS-3 altimeter data with independent measurements 
of significant waveheight made by NOAA Data Buoy Office buoys. The 
algorithm appears to be valid for a variety of ocean conditions, 
since the GEOS-3 significant waveheight estimates and the buoy sig- 
nificant waveheight measurements agree quite well for both low and 
high sea states. 

Analysis of the altimeter estimates and the buoy measure- 
ments of significant waveheight indicates that a mean bias of about 
41 cm exists between the two determinations with the GEOS-3 esti- 
mates being higher. This is consistent with a 15iS buoy overcorrec- 
tion for noise. The bias could easily be elir'"*nated by adjustment 
of the calm sea risetlme parameter in the altimeter significant 
waveheight algorithm. However, because such a small difference is 
insignificant for all known applications and because a large amount 
of GEOS-3 data has been distributed with the calm sea risetime coef- 


U)(i 



fiiMrtit ufUHl If; this invo.st » fiot ion , no chnngfrs were made to the 
algor ithm used in the teduclion and distribution (*f the data* 

Linear r*'^gfession analysis performed upon the GE0S~3 sig- 
nificant waveheight estimates and the NOAA buoy significant 
waveheight measurements indicate a high degree of correlation. 
Furthermore* the standard error of estimate was essontially 
unchanged when the 0K0S--3 estimate was changed from the dependent 
variable to the independent variable, 'nils analysis indicates that 
the relative accuracies of the altimeter estimates and the buoy 
measurements of significant waveheight are about the same. 

life results presented in Chapter ^ clearly demonstrate the 
accuracy and the reUablllty of the SWM estimation algorithm derived 
in this investigation. Additionally , the design specification of + 
20% accuracy for 2.0 < SWH < 10 meters was met. Furthermore, 
because estimates can be computed in near-real-time, sigrificant 
waveheight measurements produced from satellite altimetry could have 
significant impact on ship routing, search and rescue operation, 
meteorological research, recreational activities and numerous other 
areas where quick and reliable sea state Information is desired. 

Although the correlation between the two determinations of 
significant waveheight has been demonstrated, additional investiga- 
tion could lead to more specific information about the relationship 
of the QKOS-3 estimates to the buoy measurements. Since the origi- 


rial set of buoy data ordered from the NOAA Data Buoy Office was 
based upon a 1/4 degree search area around the buoy and since it has 
been demonstrated that the search area could be made larger than one 
degree, the number of comparisons could be greatly enlarged from 125 
if the buoy data search area was enlarged. 

If a much larger set of buoy data was available, the rela- 
tive accuracy of the altimeter estimates to the buoy measurements of 
significant waveheight could be studied more extensively. In par- 
ticular, the buoy data could be separated by buoy type. In this 
way, not only could a more powerful statement be made about the 
accuracy of the GEOS-3 significant weveheight estimates, but the 
relative accuracies of the different types of buoys could be more 
accurately determined. 

A data set consisting of 17 ground track crossings of 
GEOS-3 and Seasat was analyzed in an attempt to verify the GOASEX 
conclusion that there existed a mean bias between the two determina- 
tions of 59 cm. That analysis concluded that the mean bias was 24 
cm. 


A global atlas of significant waveheight contours was 
presented. This atlas, which contained 1500 hours of GEOS-3 data, 
took several months of computer time to compile. Contour maps were 
presented for the months December to February, March to May, June to 
August and September to November. Additionally, contour maps for 
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tho ('omposite of the entire GKOS-3 mi.ssion wore presented. 

The GEOS-3 contour maps were compared with the contour 
maps in the U. S. Navy Climatic Atlas and significant differences 
were discovered. That these two atlases differed substantially la 
an important result in itself and indicates the need for future 
study. In general, the GEOS-3 estimates of significant waveheight 
showed a higher percentage of high sea states and a lower percentage 
of low Sea states then did the Navy estimates. 

The obvious practicality of significant waveheight esti- 
mates made on a global, near-real-time basis from space is currently 
offset by the facts that only GEOS-3 and Seasat were equipped to 
make such estimates and that, due to onboard power constraints, the 
GEOS-3 altimeter was turned on and off at scheduled times. Since 
the on/off schedule was determined several weeks in advance, it was 
often difficult to obtain significant waveheight estimates at the 
time and place they were needed. Additionally, it is evident that a 
single satellite with altimeter measurements only at nadir will not 
supply the capability to produce real-time estimates of significant 
Waveheight with global coverage. These shortcomings are yet further 
arguments for a series of satellites, equipped with altimeters which 
operate continuously. Such a series of satellites would provide 
global estimates of significant waveheight bo the scientific, mili- 
tary and industrial communities. 
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APPENDIX A 


GLOBAL ATLAS OF GEOS-3 SIGNIFICANT WAVEHEIGHT DATA 


FlGUWi PAGE 

A.l ftjcember Hirough February - Percent Significant 

Wavehelght <1,5 Meters ...... . , l 

A. 2 Deceml)er nirough February > Percent Significant 

Waveheight <2.5 Meters 2 

A. 3 December Through February - Percent Significant 

Waveheight > 3.5 Meters . 3 

A. 4 December llirough February * Percent Significant 

Waveheight >6.0 Meters 4 

A. 5 ^iarch TJvrough May - Percent Significant 

Waveheight <1.5 Nteters . . 5 

A. 6 Mardfi Tlirough May - Percent Significant 

Waveheight <2.5 Meters ... 6 

A. 7 March Tlvrough May - Percent Significant 

Waveheight > 3.5 Meters 7 

A. 8 March llirough May - Percent Significant 

Waveheight > 6.0 Meters ............ 8 

A. 9 vJune Tlirough August - Percent Significant 

Waveheight <1.5 Met'^rs 9 

A. 10 June Tlirough August - Percent Significant 

Waveheight <2.5 Meters 10 

A. 11 June Tlirough August - Percent Significant 

Waveheight > 3.5 Meters 11 

A, 12 June Tlirough August - Percent Significant 

Waveheight > 6.0 Meters 12 

A. 13 September 'Flirough November - Percent Significant 

Waveheight < 1.5 Meters .... ....... 13 

A. 14 September Through November - Percent Significant 

Waveheight < 2.5 Meters 14 

A. 15 September Tlirough November - Percent Significant 

Waveheight > 3.5 Meters 15 

A. 16 September llirough November - Percent Significant 

Waveheight > 6.0 Meters 16 



FIGURE PAGE 

A. 17 Entire Mission - Percent Significant 

Waveheight < 1.5 Meters .... 17 

A. 18 Entire Mission » Percent Significant 

Waveheight <2.5 Meters 18 

A, 19 Entire Mission - Percent Significant 

Waveheight > 3.5 Meters 19 

A. 20 Entire Missiwi - Percent Significant 

Waveheight > 6.0 Meters ..... , 20 
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state contour maps are presented, and the data are displayed both by seasons and for the entire 
duration of the GEOS-3 mission. The contour maps are then compared with contour maps compiled 
by the U. S. Navy and significant differences are found. 
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